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Abstract. 

Research  on  interactions  between  quantum  gases  and  impurities  (ions,  other  atoms  etc.)  has  become  an  active  area  of 
AMO  research  over  the  past  couple  of  years.  Among  the  topics  of  interest  one  finds  the  formation  of  mesoscopic 
ions,  charge  hopping  and  quantum  transport,  self  trapping  of  impurities  in  Bose-Einstein  Condensate  (BECs),  and 
the  effect  of  impurity  particles  on  the  BEC  density  distribution.  To  study  ion-BEC  interactions,  a  method  must  be 
used  that  allows  one  to  hold  a  cold  ion  in  a  BEC  for  a  sufficiently  long  time  that  these  phenomena  become 
observable,  in  a  way  that  does  not  heat  the  ion. 

In  the  project  two  setups  have  been  constructed  and  operated  towards  these  goals.  In  the  first,  more  basic  one,  a 
metal  tip  (diameter  125  microns)  is  mounted  close  to  a  magneto-optic  trap  (MOT).  Rydberg  excitations  are 
implanted  in  a  cold-atom  region  (diameter  about  10  microns,  length  about  50  microns).  Interactions  between  the 
excitations  cause  excitation  blockades,  short-  and  medium-range  order,  as  well  as  dynamics  due  to  van-der-Waals 
interaction.  All  these  topics  have  been  studied.  Two  papers  were  published  [one  PRL,  one  PRA,  one  submitted  to 
PRL]. 

The  second  setup  is  a  double-trap  system,  in  which  we  prepare  magnetically  trapped  rubidium  atom  samples 
between  1  and  10  micro-Kelvin  and  densities  of  about  10A12  cmA-3.  We  have  prepared  and  trapped,  for  the  first 
time,  cold  circular  Rydberg  atoms  in  these  traps.  In  addition  to  being  able  to  study  the  dynamics  of  these  highly 
dipolar  atoms  in  the  magnetic  trap,  we  also  observed  a  collision-induced  loss  we  have  traced  back  to  collisions 
between  the  long-lived  circular  atoms  and  cold  ground  state  atoms.  These  types  of  collisions  are  unexplored,  to  date. 
In  this  setup  we  also  have  prepared  and  studied  D-type  molecules  between  Rydberg  and  ground-state  atoms.  These 
are  due  to  an  attractive  triplet  s-wave  low-energy  scattering  interaction  between  the  Rydberg  electron  and  a  ground- 
state  atom.  We  have  published  two  papers,  one  in  PRL,  on  in  PRA. 


Archival  publications  during  reporting  period  (all  attached): 

[*  ]  “Measurement  of  van-der-Waals  interaction  by  atom  trajectory  imaging,”  N.  Thaicharoen,  A.  Schwarzkopfy, 
and  G.  Raithel  (submitted  to  Phys.  Rev.  Lett.,  June  2015). 

[1]  “Angular-momentum  couplings  in  long-range  Rb2  Rydberg  molecules,”  D.  A.  Anderson,  S.  A.  Miller,  G. 
Raithel,  Phys.  Rev.  A  90,  062518  (2014). 

[2]  “Photoassociation  of  long-range  nD  Rydberg  molecules,”  D.  A.  Anderson,  S.  A.  Miller,  and  G.  Raithel,  Phys. 
Rev.  Lett.  112,  163201  (2014). 

[3]  “Spatial  correlations  between  Rydberg  atoms  in  an  optical  dipole  trap,”  A.  Schwarzkopf,  D.  A.  Anderson,  N. 
Thaicharoen,  and  G.  Raithel,  Phys.  Rev.  A  88 , 061406(R)  (2013);  Editor's  suggestion. 

[4]  “Production  and  trapping  of  cold  circular  Rydberg  atoms,”  D.  A.  Anderson,  A.  Schwarzkopf,  R.  E.  Sapiro,  and 
G.  Raithel,  Phys.  Rev.  A  88 , 031401  (2013);  Editor's  suggestion. 

[5]  “Imaging  Spatial  Correlations  of  Rydberg  Excitations  in  Cold  Atom  Clouds,”  A.  Schwarzkopf,  R.  E.  Sapiro,  and 
G.  Raithel,  Phys.  Rev.  Lett.  107,  103001  (2011). 


Activity  (1):  We  have  imaged  spatial  correlations  of  Rydberg  excitations  in  cold  atom  clouds.  We  use  direct  spatial 
imaging  of  cold  85Rb  Rydberg  atom  clouds  to  measure  the  Rydberg-Rydberg  correlation  function.  Initially  obtained 
results  were  in  qualitative  agreement  with  theoretical  predictions  elsewhere.  We  determined  the  blockade  radius  for 
several  D-level  Rydberg  states  and  investigated  the  dependence  of  the  correlation  behavior  on  excitation  conditions 
and  detection  delay.  Experimental  data  hinted  at  the  existence  of  long-range  order.  We  then  extended  the  research  to 
include  light- shift  potentials  generated  by  an  optical  dipole  trap.  The  optical  dipole  trap  is  used  to  prepare  dense 
ground-state  atom  samples  and  to  add  light-shifts  to  the  ground-Rydberg  transition.  We  studied  and  interpreted  the 
effects  of  laser  detuning  and  light-shifts  on  the  Rydberg  pair-correlation  function.  The  result  showed  an 
enhancement  of  the  probability  of  exciting  two  Rydberg  atoms  at  a  particular  separation  distance,  which  is 
interpreted  as  a  result  of  direct  two-photon  excitation  of  Rydberg  atom  pairs.  In  another  experiment,  we  investigated 
the  effect  of  the  time  delay  between  the  excitation  pulses  and  the  read-out  sequence.  We  observed  a  time 
dependence  of  the  pair-correlation  signal  which  is  consistent  with  the  effect  of  van-der-Waals  forces  acting  between 
Rydberg-atom  pairs.  The  results  are  significant  because  of  the  wide  use  of  optically  trapped  Rydberg  atoms  in  both 
ongoing  and  in  proposed  research.  Three  papers  (one  submitted,  two  published)  are  attached  [*,3,5]. 

Activity  (2):  We  have  successfully  prepared  and  trapped  circular  Rydberg  atoms.  Circular  Rydberg  atoms  exhibit  a 
unique  combination  of  extraordinary  properties:  long  lifetimes  (~nA5),  large  magnetic  moments  and  angular 
momenta  (|m|=4=n-l),  and  no  first  order  Stark  shift.  Here,  n,  1  and  m  are  the  principal,  orbital  and  magnetic  quantum 
numbers,  respectively.  These  features  have  made  circular  Rydberg  atoms  attractive  for  a  number  of  applications.  We 
realized  a  magnetic  trap  for  laser-cooled  Rydberg  atoms  that  are  initially  in  a  well-defined  circular  Rydberg  state. 
The  Rydberg-atom  trap  was  characterized  using  state- selective  electric-field  ionization,  direct  spatial  imaging  of  the 
atom  distributions,  and  time-of-flight  analysis  of  the  ion  signal.  At  room  temperature,  we  observed  70  percent  of  the 
trapped  atoms  remaining  after  6ms.  We  measured  an  enhancement  of  the  center-of-mass  trap  oscillation  frequency 
by  the  expected  factor  of  sqrt(m).  Simulations  of  the  state-evolution  of  circular- state  atoms  in  our  magnetic  trap, 
held  at  300K  radiation  temperature,  were  performed  and  results  were  in  good  agreement  with  the  observed 
phenomena.  One  published  paper  is  attached  [4]. 

Activity  (3):  One  of  the  most  prominent  recent  observations  of  impurities  in  dense  cold  atom  clouds  is  that  of 
Rydberg  molecules.  A  recently  discovered  class  of  long-range  Rydberg  molecules  generated  through  an  attractive 
Fermi  potential  has  generated  a  significant  amount  of  theoretical  and  experimental  interest.  Here,  we  have  observed 
long-range  homonuclear  diatomic  nD  Rydberg  molecules  photoassociated  out  of  an  ultracold  gas  of  87Rb  atoms  for 
34<n<40.  We  measured  the  ground-state  binding  energies  of  87Rb  (nD+5S_l/2)  molecular  states  to  be  larger  than 
those  of  their  87Rb(nS+5S_l/2)  counterparts,  showing  the  dependence  of  the  molecular  bond  on  the  angular 
momentum  of  the  Rydberg  atom.  We  probed  the  transition  of  87Rb(nD  +  5S_l/2)  molecules  from  the  molecular- 
binding-dominant  regime  at  low  n  [Hund’s  cases  (a)]  to  a  fme-structure-dominant  regime  at  high  n  [Hund's  case  (c)]. 
A  Fermi  model  that  includes  the  fine  structure  of  the  nD  Rydberg  atom  and  hyper  fine  structure  of  the  5S1/2 
pertuber  is  presented  that  describes  this  transition.  In  ongoing  work,  we  include  3S,  3P,  IS  and  IP  scattering 
channels  for  accurate  modeling  of  the  molecular  spectra.  Two  published  papers  are  attached  [1,2]. 


Measurement  of  van-der-Waals  interaction  by  atom  trajectory  imaging 


N.  Thaicharoen*,  A.  Schwarzkopf'' .  and  G.  Raithel 
Department  of  Physics,  University  of  Michigan,  Ann  Arbor,  MI  48109 
(Dated:  June  6,  2015) 

We  study  the  repulsive  van  der  Waals  interaction  of  70*Si/2  cold  rubidium  Rydberg  atoms  by 
analysis  of  time-delayed  pair  correlation  functions.  After  excitation,  Rydberg  atoms  are  allowed  to 
accelerate  under  the  influence  of  the  van  der  Waals  force,  before  measuring  their  positions  using  a 
single- atom  imaging  technique.  From  the  average  pair  correlation  function  of  the  atom  positions 
we  obtain  the  initial  atom-pair  separation  and  the  terminal  velocity,  which  yield  the  van  der  Waals 
interaction  coefficient  Cq.  The  measured  Cq  value  agrees  well  with  calculations.  The  experimental 
method  has  been  validated  by  simulations.  The  data  hint  at  anisotropy  in  the  overall  expansion, 
caused  by  the  shape  of  the  excitation  volume.  Our  measurement  implies  that  the  interacting  entities 
are  individual  Rydberg  atoms,  not  groups  of  atoms  that  coherently  share  a  Rydberg  excitation. 

PACS  numbers:  32.80.Ee,  34.20.Cf 


The  van  der  Waals  interaction  is  important  in  the  de¬ 
scription  and  control  of  interactions  in  few-  and  many- 
body  dynamics  studies.  This  interaction  has  been  crit¬ 
ical  in  the  observation  of  Rydberg  excitation  blockades 
and  collective  excitations  [1-4],  Rydberg  crystals  [5,  6], 
and  Rydberg  aggregates  [7,  8].  Rydberg  interactions 
have  been  used  in  quantum  information  processsing  [9- 
12].  The  van  der  Waals  interaction  between  two  Ryd¬ 
berg  atoms  has  been  measured  using  spectroscopic  meth¬ 
ods  [13,  14].  Several  measurements  have  been  per¬ 
formed  near  surfaces  to  observe  radiative  Rydberg-level 
shifts  caused  by  image  charge  interaction  near  metal  sur¬ 
faces  [15,  16].  The  van  der  Waals  interaction  between 
excited  cesium  atoms  and  a  dielectric  surface  has  been 
measured  using  selective  reflection  spectroscopy  [17]. 

Here,  we  develop  a  method  to  study  the  van  der  Waals 
interaction  between  Rydberg  atoms  using  direct  spatial 
imaging  of  their  trajectories  [18-21].  Pairs  of  S i/2  ru¬ 
bidium  Rydberg  atoms  are  prepared  with  a  well-defined 
initial  separation  by  detuning  an  excitation  laser  and  uti¬ 
lizing  the  r~ 6  dependence  of  the  van  der  Waals  interac¬ 
tion  [21,  22].  After  preparation,  the  atoms  are  subject  to 
van  der  Waals  forces  (which  are  repulsive  in  this  case). 
The  effect  of  the  forces  is  observed  by  tracking  the  inter¬ 
atomic  distance  between  the  Rydberg  atoms,  after  they 
have  been  allowed  to  move  for  selected  wait  times  (see 
Fig.  1).  The  atom  trajectories  and  thereby  the  van  der 
Waals  interaction  coefficient  Cq  are  extracted  from  the 
pair  correlation  functions  of  the  Rydberg  atom  positions. 

The  experimental  setup  is  shown  in  Fig.  1(a).  85  Rb 
ground-state  atoms  are  prepared  in  a  magneto-optical 
trap  (MOT)  at  a  density  of  >  1010  cm-3.  The  two- 
photon  Rydberg  excitation  to  70Si/2  is  driven  by  simul¬ 
taneous  780  nm  and  480  nm  laser  pulses  with  a  5  /is 
duration  and  ^1  GHz  red-detuning  from  the  5P3/2  inter¬ 
mediate  state.  Both  beams  propagate  in  the  xy  plane  and 
are  linearly  polarized  along  z.  The  780  nm  beam  has  a 
Gaussian  beam  parameter  wq  of  0.75  mm  and  the  480  nm 
beam  is  focused  to  wq  =  8  /am.  The  Rydberg  atoms  are 
ionized  by  applying  a  high  voltage  to  a  tip  imaging  probe 
(TIP)  with  a  rounded  tip  of  diameter  125  /am.  Ions  are 


accelerated  by  the  TIP  electric  field  towards  a  microchan- 
nel  plate  (MCP).  Ion  impacts  result  in  blips  produced 
by  the  MCP-phosphor  detector  assembly,  revealing  the 
Rydberg  atom  positions.  In  each  experimental  cycle  we 
record  one  image,  taken  by  a  CCD  camera,  that  typically 
contains  several  blips.  For  more  detail  see  reference  [21]. 

The  excitation  volume  is  ~470  /am  above  the  tip, 
which,  in  combination  with  the  radial  divergence  of  the 
TIP  electric  field  and  the  MCP  front  voltage,  results  in 
a  measured  magnification  of  155  with  an  uncertainty  of 
2%.  We  have  calibrated  the  magnification  by  project¬ 
ing  a  known  double-slit  interference  pattern  of  480  nm 
excitation  light  into  the  object  plane  and  measuring  the 
distribution  of  counts  in  the  image  (MCP)  plane.  The 
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FIG.  1.  Geometry  (a)  and  timing  sequence  (b)  of  the  exper¬ 
iment.  The  wait  time  is  the  time  difference  between  the  end 
of  the  excitation  and  the  onset  of  the  field  ionization  pulse, 
(c)  Van  der  Waals  potential  versus  separation  between  two 
Rydberg  atoms.  The  separation  increases  due  to  the  repul¬ 
sive  van  der  Waals  force.  We  image  the  resulting  trajectories 
by  varying  the  wait  time.  In  the  insets  in  (c),  circles  represent 
two  Rydberg  atoms  evolving  under  the  influence  of  the  van 
der  Waals  force.  The  kinetic  energy  release,  indicated  by  the 
blue  velocity  arrows,  equals  the  van  der  Waals  potential  Wo 
at  the  initial  separation  ro. 
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magnification  uncertainty  follows  from  the  variation  of 
the  peak  positions  from  an  equidistant  series  (which  one 
would  expect  for  the  double-slit).  The  resolution  in  the 
object  plane  is  ~  1  /i m;  it  follows  from  the  magnification 
and  the  point  spread  function  of  the  detector.  We  can 
easily  resolve  the  correlation  length  between  the  Rydberg 
atoms,  which  is  on  the  order  of  10  /i m.  The  field  of  view 
typically  contains  5  to  10  detected  atoms.  In  the  image 
processing  sequence,  we  first  extract  ion  impact  positions 
on  the  MCP  by  using  a  peak  detection  algorithm  [21]. 
Out  of  10000  images  taken  in  each  dataset,  we  select  the 
5000  images  with  the  highest  numbers  of  detected  ions. 
From  this  subset,  we  calculate  the  average  pair  correla¬ 
tion  image  and  normalize  it  such  that  at  large  distances 
it  approaches  the  value  of  one. 


2.5|is  IOjis  20(is  30(is  4(+s 


FIG.  2.  Pair  correlation  images  from  experiment  (top  row) 
and  simulation  (bottom  row)  for  the  selected  wait  times.  The 
linear  grayscale  ranges  from  0  (white)  to  2  (black).  Values 
of  1,  <1,  and  >1  indicate  no  correlation,  anticorrelation,  and 
positive  correlation,  respectively. 

The  excitation  of  the  7051/2  Rydberg  level  is  detuned 
by  =  4  ±  2  MHz  with  respect  to  the  two-photon  reso¬ 
nance.  The  positive  excitation  detuning  is  compensated 
by  the  (repulsive)  van  der  Waals  interaction.  We  prefer¬ 
entially  excite  Rydberg  atoms  in  pairs  at  a  separation  at 
which  the  van  der  Waals  interaction  is  25ui,  =  8+4  MHz. 
We  determine  ro  from  pair-correlation  images.  After  ex¬ 
citation,  the  Rydberg  atoms  are  allowed  to  move  within 
selected  times  before  measuring  their  positions  by  apply¬ 
ing  a  field  ionization  pulse  [Fig.  1(b)].  The  atom-position 
data  are  processed  as  explained  above,  yielding  average 
pair  correlation  images.  In  Fig.  2,  we  present  the  pair 
correlation  images  for  wait  times  2.5  /i s,  10  /i s,  20  /i s, 
30  /is,  and  40  /is.  The  average  interaction  time  of  the 
Rydberg  atoms  is  the  wait  time  plus  half  the  excitation 
pulse  length  (wait  time  +  2.5  /i s).  The  pair  correlation 
image  at  2.5  /is  exhibits  strong  correlation  enhancement 
at  a  fairly  well-defined  radius.  The  initial  correlation 
is  critical  for  our  trajectory  experiment;  it  is  sufficient 
to  track  Rydberg-pair  trajectories  out  to  approximately 
40  /i s.  With  increasing  wait  time,  the  radius  of  enhanced 
pair  correlation  increases,  reflecting  an  increase  of  inter¬ 
atomic  separation  due  to  the  repulsive  van  der  Waals 
interaction.  At  long  wait  times,  the  enhancement  ring  is 
blurred  out  due  to  initial  thermal  atom  velocities  in  the 
MOT  (temperature  ~  100  /iK). 

To  determine  the  most  probable  separation  rp(£)  be¬ 


tween  Rydberg- atom  pairs  at  wait  time  £,  we  first  com¬ 
pute  the  angular  integrals  /(r )  of  the  experimental  pair 
correlation  images  in  Fig.  2,  as  shown  in  Fig.  3(a).  The 
values  of  rp(£)  are  obtained  from  local  parabolic  fits  cen¬ 
tered  approximately  at  the  peak  positions  of  the  /(r ) 
curves.  We  include  4  to  7  data  points  in  the  fit  (depend¬ 
ing  on  the  shapes  of  the  curves).  The  resulting  separa¬ 
tions  rp(£)  are  shown  in  Fig.  3(b).  The  visibility  of  the 
pair  correlation  enhancement,  also  shown  in  Fig.  3(b), 

is  (/max  /min) /(/max  +  /min) 5  where  /max  IS  the  peak 

value  of  /(r)  in  the  range  >  10  /im  and  /m in  is  the  min¬ 
imal  value  near  7  /im.  A  cursory  inspection  of  Fig.  3(b) 
already  shows  that  the  trajectory  of  the  Rydberg- atom 
pairs  is  characterized  by  an  initial  acceleration  phase, 
during  which  the  initial  van  der  Waals  potential  energy, 
Wo,  is  converted  into  kinetic  energy,  and  a  later  phase 
during  which  the  atoms  keep  separating  at  a  fixed  veloc¬ 
ity.  The  drop  in  visibility  is  due  to  the  thermally-induced 
blurring  of  the  correlation  ring  at  late  times. 


FIG.  3.  (a)  Angular  integrals  /(r)  of  the  pair  correlation 

images  in  Fig.  2  at  wait  times:  2.5  /is  (red  solid),  10  /is  (or¬ 
ange  dashed),  20  /is  (green  dotted),  30  /is  (blue  dot-dashed), 
and  40  /is  (pink  short- dashed).  The  y  axis  is  for  the  red- 
solid  curve;  for  clarity,  the  other  curves  are  shifted  down  in 
equidistant  intervals  of  0.3.  (b)  Interatomic  separations  be¬ 
tween  Rydberg-atom  pairs  as  a  function  of  interaction  time 
(left  axis,  blue  circles),  obtained  from  the  peak  positions  in 
the  /(r)  curves.  The  blue  dashed  line  represents  a  linear  fit 
at  long  wait  times  (20-40  /is).  The  red  solid  line  shows  sim¬ 
ulation  results  obtained  for  =  3  MHz.  The  pink  hollow 
squares  show  the  visibility  of  the  experimental  pair  correlation 
enhancement  (right  axis)  as  defined  in  text. 

To  extract  the  van  der  Waals  Cq  coefficient,  one  may 
consider  an  isolated  atom  pair  excited  at  an  initial  sepa¬ 
ration  ro.  The  initial  van  der  Waals  energy  is 

Wo  =  %  =  2  6ul.  (1) 

ro 

From  Eq.  1,(76  =  2£z/l^o  carb  in  principle,  be  obtained 
from  spectroscopic  measurement  [14].  This  method  re¬ 
quires  a  well-defined  ro,  a  narrow  laser  linewidth,  and 
accurate  knowledge  of  Su l-  In  our  work,  the  relative  un¬ 
certainty  in  rp  is  about  12%  (because  the  relative  mag¬ 
nification  uncertainty  is  2%).  In  comparison,  the  uncer¬ 
tainty  arising  from  =4  +  2  MHz  is  much  larger.  Tak¬ 
ing  all  uncertainties  into  account,  Eq.  1  leads  to  values  of 
Cq  ranging  from  7xl0-58  Jm6  to  40xl0-58  Jm6.  There- 
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fore,  Eq.  1  only  allows  us  to  perform  a  crude  order-of- 
magnitude  estimate  for  Co  .  The  main  weakness  of  Eq.  1 
is  that  the  C^-values  derived  from  it  directly  reflect  our 
large  relative  uncertainty  in  Jz/l- 

A  better  way  for  us  to  obtain  Cq  is  to  track  the  evolu¬ 
tion  of  the  Rydberg-atom  trajectories.  Over  the  exper¬ 
imentally  investigated  wait  times,  the  entire  initial  van 
der  Waals  energy  Wo  becomes  converted  into  kinetic  en¬ 
ergy.  With  the  reduced  mass  of  the  atom  pair,  /q  and 
the  terminal  relative  velocity,  vt,  it  is 

Wo  =  ^  =  (2) 

'0  z 

To  obtain  ro,  we  take  the  weighted  average  of  the  fit  re¬ 
sults  rp(t)  at  the  earliest  times  used  in  the  experiment 
[the  points  in  Fig.  3(b)  at  250  ns  and  2.5  y s].  The  aver¬ 
aging  is  valid  because  during  the  first  few  microseconds 
after  excitation  the  Rydberg  atoms  are  frozen  in  place 
due  to  their  inertia.  The  statistical  weights  are  given  by 
the  inverse  squares  of  the  fitting  uncertainties  of  rp  at 
250  ns  and  at  2.5  ys.  We  obtain  ro  =  8.78  y m  with  a  net 
fitting  uncertainty  of  0.02  /im.  To  determine  t?t,  we  per¬ 
form  a  linear  fit  at  long  wait  times  (20-40  ys)  and  obtain 
vt  =  0.182  m/s  with  a  fitting  uncertainty  of  0.008  m/s.  It 
follows  Cq  =  |/rr2 Tq  =  (5.4T0.5)  x  10-58  Jm6.  Including 
the  2%  magnification  uncertainty,  the  total  uncertainty 
of  Cq  becomes  1.0  xlO-58  Jm6. 

The  final  relative  uncertainty  of  Co  follows  from  three 
statistically  independent  contributions:  the  magnifica¬ 
tion  uncertainty,  the  fit  uncertainty  for  rt,  and  the  fit 
uncertainty  for  ro-  The  respective  powers  at  which  these 
quantities  enter  into  Co  are  8,  2,  and  6.  Factoring  in 
these  powers,  the  three  quantities  contribute  respective 
independent  uncertainties  of  16%,  8.7%,  and  1.4%  to  the 
relative  uncertainty  of  Co ,  leading  to  the  total  relative 
uncertainty  of  18%.  In  our  method,  achieving  a  small 
magnification  uncertainty  is  particularly  important. 

A  semi-classical  3D  simulation  of  the  dynamics  of  Ryd¬ 
berg  atoms  interacting  by  an  isotropic  van  der  Waals 
force  has  been  performed  to  interpret  the  experiment. 
The  simulation  volume  of  (140  /im)3  exceeds  the  exper¬ 
imentally  relevant  volume  by  about  a  factor  of  two  in 
each  dimension.  The  480  nm  excitation  beam  propagates 
along  the  y  direction.  The  excitation  volume  in  the  trans¬ 
verse  directions  ( x  and  z)  is  limited  by  the  size  of  the  ex¬ 
citation  beam  (wo  =  8  /im).  The  number  of  simulations 
is  5000,  in  keeping  with  the  number  of  images  analyzed  in 
the  experiment.  Rydberg-atom  positions  and  velocities 
are  initialized  as  explained  below,  and  the  positions  and 
velocities  are  then  propagated  using  a  Runge-Kutta  inte¬ 
grator  that  includes  all  pair-wise  interatomic  forces.  We 
use  Co  =  5.7  x  10-58Jm6  [23]  in  the  initialization  and 
integration  procedures.  To  avoid  edge  effects,  we  crop 
the  xy  processing  area  in  the  images  to  (70  /im)2  before 
calculating  the  average  pair-correlation  images  and  the 
radial  functions  /(r). 

In  each  simulation,  we  start  by  drawing  the  number 
of  Rydberg  atoms  to  be  used,  7iRyd,  from  a  Poissonian 


distribution.  Initial  trial  positions  of  the  Rydberg  atoms 
are  then  drawn  from  a  probability  distribution  that  is 
Gaussian  along  x  and  z,  with  a  wo  of  8  /im,  and  uniform 
along  y  (in  close  analogy  with  the  experiment).  For  an 
atom  i  at  a  trial  position  iq,  the  effective  detuning  Jeff  (ri) 
due  to  all  other  atoms  j  at  positions  Yj  which  are  already 
excited  into  the  Rydberg  state  is 

<5eff(ri)  =  y]  |r.  _  r_.|6  -  ^  (3) 

Since  the  detuning  is  substantial,  the  first  pair  of 
Rydberg  atoms  is  simultaneously  excited  via  off-resonant 
excitation  [21]  (because  the  intermediate  state  in  which 
there  is  only  one  Rydberg  atom  present  is  off-resonant). 
Therefore  the  value  of  S  in  Eq.  3  for  the  first  Rydberg- 
atom  pair  is  set  to  be  twice  the  laser  detuning  S z/l-  For 
the  excitation  of  additional  atoms  we  use  S  =  Jz/l  to 
simulate  the  stepwise  addition  of  those  atoms,  which 
can  be  a  near-resonant  process  with  other  atoms  already 
present  [8] .  Detunings  due  to  the  Doppler  effect  are  about 
300  kHz  and  are  neglected.  The  excitation  probabil¬ 
ity  Pex(^eff)  is  a  Gaussian  centered  at  5e r  =  0  with  a 
FWHM  of  4  MHz  (given  by  the  excitation  bandwidth) 
and  Pex(^eff  =  0)  =  1.  A  new  Rydberg  atom  at  position 

is  created  if  the  excitation  probability  is  larger  than  a 
number  randomly  drawn  between  0  and  1.  This  proce¬ 
dure  is  repeated  with  new  trial  positions  until  the  desired 
number  of  7iRyd  Rydberg  atoms  has  been  reached.  The 
initial  center-of-mass  velocities  of  the  atoms  are  assigned 
using  a  Maxwell  distribution  at  temperature  100  yK. 

In  the  simulation,  we  record  the  atom  positions  at  the 
same  wait  times  as  used  in  the  experiment.  The  resulting 
pair  correlation  images  are  shown  for  the  case  Jz/l  = 
3  MHz  in  the  bottom  row  in  Fig.  2.  Black  pixels  along  the 
left  and  right  edges  of  each  pair  correlation  are  an  artifact 
due  to  the  normalization  used  in  the  image  processing. 

From  the  simulated  pair  correlation  images  we  cal¬ 
culate  the  I(r )  curves,  as  in  the  experiment,  in  order 
to  obtain  the  most  probable  separations  rp(£)  between 
Rydberg-atom  pairs  at  each  interaction  time  [see  solid 
curve  in  Fig.  3(b),  which  is  for  =  3  MHz].  The 
simulated  and  experimental  results  for  rp(£)  are  in  good 
agreement.  In  order  to  test  how  well  our  experimental 
procedure  reproduces  the  Co  coefficient  that  underlies 
the  atomic  interactions,  we  evaluate  the  simulated  results 
for  rp(£)  using  the  same  method  that  we  also  employ  for 
the  experimental  data.  The  results  for  ro,  vt  and  Cq 
extracted  from  the  simulated  rp(£)  curves  are  shown  in 
Table  I  for  four  choices  of  S These  Co  values  agree, 
within  the  uncertainties,  with  the  value  that  has  been 
entered  as  a  fixed  input  into  the  simulation.  This  finding 
validates  the  experimentally  used  procedure;  in  particu¬ 
lar  it  is  seen  that  the  method  is  not  very  sensitive  on 
The  experimental  and  calculated  [23]  Co  values  are  also 
included  in  Table  I  for  reference. 

The  experimental,  simulated,  and  expected  values  for 
Co  in  Table  I  are  in  reasonable  agreement.  Moreover,  we 
observe  that  the  entire  simulated  curve  rp(£)  for  Jz/l  = 
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TABLE  I.  Comparison  of  calculated,  experimental,  and  sim¬ 
ulated  results  for  Cq,  7*0,  and  vt. 


ro 

W 

Vt 

(m/s) 

c6 

(xlO-58Jm6) 

Calculation 

Experiment 

8.8T0.2 

0.18T0.09 

5.7 

5.4T1.0 

Simulations 

=  2  MHz 

8.85T0.06 

0.17T0.02 

4.9  ±1.0 

5is l  —  3  MHz 

8.14T0.03 

0.233±0.005 

5.5±0.3 

—  4  MHz 

7.76T0.05 

0.278T0.003 

6.0±0.3 

=  5  MHz 

7.51T0.04 

0.28T0.01 

4.9±0.5 

FIG.  4.  Pair  correlation  function  at  wait  time  30  ys  from 
(a)  experiment  and  (b)  simulation  at  =  3  MHz.  The 
anisotropic  expansion  of  the  atom  ensemble  causes  the  radius 
of  enhanced  pair  correlation  along  the  transverse  direction  (T) 
to  be  larger  than  along  the  excitation  beam  direction  (y). 


3  MHz  matches  the  experimental  result  very  well  [see 
Figs.  2  and  3(b)].  Overall,  our  findings  lend  credibility 
to  our  method  of  measuring  Cq. 

We  note  that  higher-order  quadrupole-dipole  and 
quadrupole-quadrupole  interaction  potentials,  which 
scale  as  r~7  and  r-8,  are  not  important  at  the  distances 
relevant  in  our  work,  at  the  current  level  of  precision.  We 
have  verified  this  in  a  calculation  of  interaction  potentials 
in  which  all  terms  up  to  the  quadrupole-quadrupole  in¬ 
teraction  have  been  included  [24]. 

A  question  of  interest  is  whether  the  interatomic  force 
is  acting  on  “superatoms”,  i.e.  Rydberg  excitations 
shared  among  a  number  of  ground-state  atoms,  or  on 
individual  atoms.  In  our  case,  there  are  several  tens  of 
atoms  within  one  blockade  sphere.  In  our  analysis,  we 
have  assumed  that  the  effective  mass  of  the  interacting 
entities  is  half  the  rubidium  atom  mass,  i.e.  we  have  im¬ 
plied  that  the  interacting  entities  are  individual  atoms 
and  not  superatoms.  For  our  density  and  blockade  ra¬ 
dius,  the  total  mass  of  a  superatom  equals  that  of  several 
tens  of  atoms.  The  agreement  between  the  Cq  coefficients 
in  Table  I  demonstrates  that  the  interacting  entities  are 
indeed  individual  atoms.  The  finding  implies  that,  dur¬ 
ing  the  course  of  the  van  der  Waals  interaction,  excita¬ 
tions  within  superatoms  become  projected  onto  individ¬ 
ual  atoms,  which  are  then  ejected  from  the  initial  super¬ 
atom  volumes.  The  phenomenon  has  been  predicted  in 
[25]  for  superatom  clouds  interacting  via  a  dipole-dipole 
interaction.  In  our  work,  we  arrive  at  a  similar  conclu¬ 
sion  for  van-der-Waals-interacting  Rydberg  atoms  in  a 
3D  system. 

Close  inspection  of  the  pair  correlation  functions  from 
the  simulation  presented  in  Fig.  2  reveals  anisotropic  ex¬ 
pansion  behavior:  the  radius  of  enhanced  pair  correlation 
is  larger  in  x  than  in  y  direction.  In  several  experimental 
data  sets  we  see  some  evidence  of  anisotropic  expansion, 
such  as  in  the  experimental  result  shown  in  Fig.  4(a). 
The  anisotropic  expansion  is  due  to  the  cylindrical  shape 
of  the  excitation  volume  (not  the  interatomic  interaction, 
which  is  isotropic).  Since  the  blockade  radius  is  close  to 
the  transverse  size  of  the  excitation  region,  at  most  two 
atoms  can  be  excited  side-by-side  in  x  direction,  lead¬ 
ing  to  free,  unimpeded  expansion  along  that  direction. 


In  contrast,  more  than  two  Rydberg  atoms  can  be  cre¬ 
ated  along  the  y  direction.  Therefore,  along  y  the  ex¬ 
pansion  is  slowed  down  due  to  multiple- atom  repulsion. 
In  our  experiment,  the  overall  Rydberg-atom  density  is 
high  enough  to  sometimes  observe  this  effect. 

In  summary,  we  have  studied  the  trajectory  of 
Rydberg-atom  pairs  interacting  by  repulsive,  isotropic 
van  der  Waals  interaction.  We  have  extracted  the  Cq 
coefficient  from  the  experiment  and  compared  it  with 
simulations  and  calculations,  and  have  observed  good 
agreement.  The  result  implies  that  the  interaction  oc¬ 
curs  between  individual  atoms,  not  superatoms.  We  have 
observed  indications  of  an  anisotropic  effect  in  the  ex¬ 
pansion,  caused  by  the  excitation  geometry.  Future  work 
may  involve  atom- atom  interactions  that  are  anisotropic, 
such  as  the  dipole-dipole  interaction. 
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Measurement  of  van-der-Waals  interaction  by  atom  trajectory  imaging 
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Supplement 

We  have  used  the  model  presented  in  detail  in  [1]  to  calculate  the  interaction  potential  between  a  pair 
of  70Si/2  rubidium  Rydberg  atoms  with  dipole-dipole  interaction  terms  only,  as  well  as  with  dipole- 
dipole,  dipole-quadrupole  and  quadrupole-quadrupole  terms.  The  projection  of  the  electron  angular 
momenta  onto  the  internuclear  axis,  M,  can  take  the  values  0  or  +/-1.  For  M=0  there  are  two  Rydberg- 
pair  70Si/2  -70Si/2  levels,  labeled  A  and  B,  while  for  M=+/-l  there  is  only  one  (see  table).  QQ  means  all 
interactions  up  to  and  including  quadrupole-quadrupole  terms  are  included.  DD  means  only  dipole- 
dipole  terms  are  included. 

It  is  seen  that  all  the  atom  pair  has  very  similar  shifts  on  all  potentials,  and  that  interaction  terms  beyond 
the  dipole-dipole  term  are  not  important  in  the  atom  distance  range  of  interest  (R>7p.m).  The  results 
also  agree  well  with  a  perturbative  level-shift  calculation  in  [2], 


M=0  Level  A 

M=0  Level  B 

M=1 

R 

QQ 

DD 

QQ 

DD 

QQ 

DD 

(Lim 

MHz 

MHz 

MHz 

MHz 

MHz 

MHz 

6 

18.215 

18.223 

19.046 

19.052 

18.424 

18.431 

7 

7.263 

7.268 

7.604 

7.609 

7.349 

7.353 

8 

3.267 

3.269 

3.421 

3.424 

3.305 

3.307 

9 

1.613 

1.614 

1.690 

1.691 

1.632 

1.633 

10 

0.858 

0.858 

0.898 

0.899 

0.868 

0.868 
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Spatial  correlations  between  Rydberg  atoms  in  an  optical  dipole  trap 
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We  use  direct  spatial  ion  imaging  of  cold  85 Rb  Rydberg  atom  clouds  to  measure  the  Rydberg-Rydberg 
correlation  function,  with  and  without  light-shift  potentials  generated  by  an  optical  dipole  trap.  We  find  that 
the  blockade  radius  depends  on  laser  detunings  and  spatially  varying  light  shifts.  At  certain  laser  detunings 
the  probability  of  exciting  Rydberg  atoms  at  particular  separations  is  enhanced,  which  we  interpret  to  be  a 
result  of  direct  two-photon  excitation  of  Rydberg  atom  pairs.  The  results  are  in  accordance  with  predictions 
[F.  Robicheaux  and  J.  V.  Hernandez,  Phys.  Rev.  A  72,  063403  (2005)]  and  a  model  we  develop  that  accounts  for 
a  one-dimensional  dipole-trap  potential. 

DOI:  10.1 103/PhysRevA.88.061406  PACS  number(s):  32.80.Ee,  32.80.Rm,  34.20.Cf,  42.50.Hz 


Rydberg-Rydberg  interactions  are  essential  in  recently 
proposed  [1,2]  and  demonstrated  [3,4]  neutral-atom  quantum 
information  processing  schemes.  They  also  play  important 
roles  in  many-body  Rydberg-atom  physics  [5-9]  and  specifi¬ 
cally  in  proposals  to  produce  Rydberg  crystals  [10-13].  Much 
recent  work  uses  or  assumes  the  presence  of  an  optical  trap 
for  ground  or  Rydberg  atoms.  Indeed,  in  many  quantum 
information  schemes  a  trap  is  important  to  assist  single-qubit 
addressability  [14-16].  Nonetheless,  to  our  knowledge  the 
influence  of  the  trap  on  the  Rydberg  blockade,  a  signature  of  the 
Rydberg-Rydberg  interactions  [17,18],  has  not  been  explicitly 
studied.  In  this  work,  we  demonstrate  that  laser  detunings 
and  background  optical  potentials  each  have  significant  effects 
on  the  Rydberg  blockade.  We  focus  on  the  blockade-induced 
correlations  between  Rydberg  atom  positions. 

We  start  with  a  magneto-optic  trap  (MOT)  of  85  Rb  atoms 
(density  >1010  cm-3)  and,  after  switching  off  the  MOT  light, 
produce  Rydberg  atoms  via  a  two-photon  excitation  from  the 
5S\/2  ground  state  by  using  780-  and  480-nm  laser  pulses.  (The 
MOT  magnetic  field  is  always  on.)  The  excitation  is  detuned 
from  the  intermediate  5P3/2  state  by  1  GHz.  We  use  crossed 
excitation  beams,  as  shown  in  Fig.  1(a).  The  780-nm  beam 
has  a  Gaussian  beam  parameter  of  wo  =  0.75  mm  while  the 
480-nm  beam  is  focused  to  /zm.  Beam  powers  are  2 

and  36  mW,  respectively.  Each  is  linearly  polarized  along  z, 
giving  Rydberg  atoms  with  mj  =  ±1/2.  In  part  of  the  work 
we  excite  the  Rydberg  atoms  in  an  optical  dipole  trap  (ODT) 
with  4  W  of  power  and  a  focal  spot  of  wo  =  23  /z m,  crossed 
with  the  480-nm  beam  focus.  The  ground-state  atom  density  in 
the  ODT  is  >10n  cm-3.  We  use  Stark  spectroscopy  to  reduce 
the  electric  field  to  <20  mV/cm.  This  is  sufficient  to  avoid 
effects  on  the  excitation  blockade  for  the  Rydberg  states  we 
excite  [19]. 

We  achieve  spatially  resolved  detection  of  the  Rydberg 
atoms  by  using  a  needle-shaped  beryllium-copper  field  ion¬ 
ization  electrode  with  a  rounded  tip  of  diameter  125  /zm 
[Fig.  1(b)],  referred  to  as  the  tip  imaging  probe  (TIP)  [18]. 
The  photoexcitation  region  is  ~220  /zm  above  the  TIP.  The 
ions  resulting  from  field  ionization  are  accelerated  along  the 
radially  divergent  electric  field  lines  produced  by  the  TIP.  The 
ions  are  detected  ~  15  cm  away  by  a  microchannel  plate  (MCP) 
spatially  sensitive  ion  detector  assembly.  The  magnification 
factor  (the  ratio  between  the  ion  separations  in  the  MCP 


plane  and  the  Rydberg-atom  separations  in  the  excitation 
plane)  is  ~300,  providing  submicron  resolution  (~0.5  /zm). 
This  is  comparable  to  existing  [20,21]  and  proposed  [22] 
diffraction-limited  optical  detection  schemes  and  scanning 
electron  microscope  schemes  [23]. 

To  measure  the  Rydberg-Rydberg  correlation  function,  we 
first  acquire  10  000  CCD  (charge-coupled  device)  images  of 
the  MCP  phosphor  screen,  each  image  showing  ion  impact 
positions  for  a  single  experimental  cycle  of  Rydberg  excitation 
and  ionization.  We  calculate  the  Rydberg  pair  correlation  for 
each  experimental  cycle  and  then  average  these  results.  We  use 
a  peak  detection  algorithm  to  detect  individual  ion  positions  in 
each  image;  we  then  calculate  the  pair  correlation  from  these 
discrete  positions.  Discrete  peak  detection  is  reliable  because 
the  minimum  separation  between  Rydberg  atoms  (the  blockade 
radius)  is  substantially  larger  than  the  imaging  resolution. 

We  measure  the  Rydberg-Rydberg  correlation  function  for 
excitation  to  states  595*1/2,  666*1/2,  745*1/2,  57D5/2,  64D5/2, 
and  72D5/2.  Laser  excitation  pulse  lengths  are  5  /zs  for  5  states 
and  1  /zs  for  D  states,  to  account  for  differences  in  oscillator 
strength  and  to  make  count  rates  comparable.  To  eliminate 


FIG.  1.  (Color  online)  Experimental  geometry:  (a)  A  top  view 
showing  the  crossed  780-  and  480-nm  excitation  beams  and  the 
1064-nm  trap  beam,  with  the  Ixjoq  widths  marked  for  the  480-  and 
1064-nm  beams.  The  circle  shows  an  estimate  of  the  field  of  view  of 
the  microchannel  plate  (MCP)  in  the  excitation  plane  (radius  25  /zm). 
(b)  A  shadow  image  of  the  ODT  above  the  tip  imaging  probe  (TIP), 
and  a  sketch  of  the  ion  trajectories  toward  the  MCP. 
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FIG.  2.  (a)  Excitation  spectrum  without  ODT  for  66Ai/2.  Note 
the  line  asymmetry,  with  interaction-induced  broadening  for  A  >  0. 
The  vertical  lines  show  frequencies  at  which  correlation  functions 
are  measured.  Panels  (b)  and  (c):  Pair  correlations  for  excitation 
of  66S1/2  at  the  indicated  detunings,  (d)  Level  diagrams  for  off- 
resonant  excitation  without  ODT.  The  laser  frequency  is  /,  and  the 
frequency  to  excite  one  Rydberg  atom  without  interactions  is  /0. 
Dotted  lines  indicate  energy  levels  without  interactions.  Solid  lines 
are  energy  levels  with  interactions.  The  diagram  is  for  the  case  of 
blue  detuning  from  /0  (A  >  0).  The  “two-photon  resonance”  for  the 
direct  excitation  of  two  Rydberg  atoms  occurs  when  the  interaction 
energy  Vmt/h  =  2 A. 


most  images  with  fewer  than  two  counts,  we  only  use  the 
5  000  (of  the  10  000  taken)  that  have  the  highest  count  numbers 
for  calculating  the  average  pair  correlation. 

For  each  state,  we  adjust  the  480-nm  laser  frequency  to 
excite  at  selected  frequency  points  on  the  excitation  spectra 
without  and  with  the  ODT,  as  indicated  by  vertical  lines 
in  Figs.  2(a)  and  3(e),  respectively.  Resulting  Rydberg  pair 
correlations  are  shown  in  Figs.  2(b),  2(c),  and  3(a)-3(d). 
Correlation  values  indicate  the  likelihood  of  finding  Rydberg 
atoms  with  a  given  projection  of  their  separation  vector  into 
the  xy  plane.  The  linear  grayscale  is  set  such  that  white  is 
0  and  black  is  2.  With  the  normalization  used,  a  value  of  1 
indicates  no  correlation,  i.e.,  the  equivalent  of  random  Rydberg 
positions.  Values  <1  indicate  anticorrelation  and  those  >1 
indicate  positive  correlation. 

We  first  perform  on-resonant  excitation  of  state  66S1/2 
without  ODT,  with  the  result  shown  in  Fig.  2(b).  For  Rydberg 
separations  r  <8.5  gm,  correlation  values  are  suppressed  due 
to  the  Rydberg  excitation  blockade.  The  single  white  pixel  in 
the  center  reflects  the  inability  to  resolve  two  Rydberg  atoms 
with  a  separation  less  than  the  imaging  resolution.  We  next 
detune  the  480-nm  laser  frequency  to  the  blue  of  resonance, 
exciting  at  the  point  where  the  count  rate  has  dropped  by 
half  [see  Fig.  2(a)].  The  resulting  pair  correlation  is  shown  in 
Fig.  2(c).  The  dark  ring  in  Fig.  2(c)  reflects  a  clear  enhancement 
of  the  probability  of  exciting  two  Rydberg  atoms  at  a  particular 
separation  distance,  as  predicted  in  Refs.  [7,9].  One  expects 
this  result  for  a  detuning  of  the  same  sign  as  the  van  der 
Waals  shift  Vmt  (positive  for  rubidium  S  states,  negative  for 
D  states  [19,24]).  The  preferred  separation  is  such  that  the 
laser  resonantly  excites  two  ground- state  atoms  directly  to 
two  Rydberg  atoms,  without  populating  an  intermediate  state 
that  carries  a  single  Rydberg  excitation.  An  energy  level 
diagram  for  this  “two-photon”  process  is  shown  in  Fig.  2(d). 
(Experimentally,  this  is  a  four-photon  transition  because  each 
Rydberg  transition  involves  both  a  780-  and  480-nm  photon.) 
For  S  states  we  find  that  blue  detunings  yield  these  correlation 
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FIG.  3.  (Color  online)  (a)  Correlation  function  for  resonant  exci¬ 
tation  of  59£i/2  without  ODT.  Panels  (b)-(d):  Correlation  functions 
for  excitation  in  the  ODT,  at  the  indicated  detunings,  (e)  Spectra 
for  excitation  of  59Ai/2  without  ODT  (black  dashed  curve)  and  with 
ODT  (red  [gray]  solid  curve).  Vertical  lines  indicate  the  excitation 
frequencies  used  to  acquire  the  correlation  functions  in  panels  (a)-(d). 
The  0-MHz  line  is  for  excitation  without  ODT;  for  the  other  lines  the 
ODT  is  on.  (f)  Angular  averages,  /(r),  for  the  correlation  functions 
in  panels  (a)-(d).  The  curves  correspond  to  correlation  functions  as 
follows:  black  solid  (a),  red  (gray)  dashed  (b),  green  (gray)  dotted 
(c),  and  purple  (dark  gray)  dot-dash  (d).  The  vertical  scale  is  labeled 
for  the  black  solid  curve;  a  vertical  offset  of  0. 1  is  applied  between 
subsequent  curves  to  allow  visual  discrimination.  For  r  less  than  the 
imaging  resolution,  each  curve  I(r)  drops  to  zero. 


enhancements,  while  red  detunings  do  not.  For  D5/2  states  we 
observe  the  opposite  behavior,  since  V{nt  <  0. 

Our  interpretation  accords  with  the  asymmetric  broadening 
of  the  excitation  lines,  such  as  the  S  lines  in  Figs.  2(a)  and  3(e) 
(dashed  line).  This  indicates  the  presence  of  the  two-photon 
process  shown  in  Fig.  2(d),  with  a  positive  Vmt/  h  of  several 
megahertz.  D5/2  lines  (not  shown)  are  broadened  on  the  low- 
frequency  side,  corresponding  to  a  negative  V[nt. 

In  subsequent  experiments,  we  excite  Rydberg  atoms  in 
the  ODT  at  selected  frequencies  on  the  ODT  spectrum,  as 
depicted  in  Fig.  3(e).  Pair  correlations  for  5 9 S1/2  are  shown 
in  Figs.  3(a)-3(d);  Fig.  3(a)  is  for  excitation  without  ODT, 
shown  for  comparison  with  ODT  results  in  Figs.  3(b)-3(d). 
With  increasing  detuning,  a  ring  of  enhanced  pair  correlation 
both  grows  in  visibility  and  shrinks  in  radius.  At  40-MHz 
detuning  a  second  ring  emerges.  We  believe  this  is  due  to 
exciting  a  string  of  three  equally  spaced  Rydberg  excitations, 
with  the  middle  excitation  at  the  center  of  the  ODT.  The 
second  ring  is  seen  especially  well  in  the  angular  averages 
7(r)  shown  in  Fig.  3(f).  The  error  bars  in  Fig.  3(f)  are  obtained 
by  testing  the  variability  of  the  results  using  subsets  of  the 
data.  The  uncertainty  increases  significantly  for  datasets  at 
large  laser  detunings  due  to  reduced  count  rates.  A  model  for 
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FIG.  4.  (Color  online)  Measurements  and  predictions  of  the 
blockade  radius  rb  for  several  nS\/2  andnD5/2  states,  for  zero  detuning 
without  ODT.  Error  bars  on  measurements  are  ±10%  of  the  measured 
value,  due  to  uncertainty  in  the  magnification  ratio  calibration. 


the  excitation  process,  which  explains  the  trends  in  Fig.  3,  is 
presented  below. 

We  measure  the  blockade  radius  rb  for  each  I(r )  curve  using 
the  point  closest  to  the  origin  where  7(r)  crosses  its  asymptotic 
value  for  large  r.  This  method  is  suggested  by  results  of  a 
stochastic,  nonquantum  simulation  of  the  excitation-blockade 
process,  which  accounts  for  our  geometry  including  the  z 
projection  during  detection.  The  hard- sphere  blockade  radii 
entered  as  inputs,  and  the  resulting  asymptote  crossings  of  the 
simulated  /(r)  curves,  agree  to  within  ±1  gm. 

In  Fig.  4  we  compare  measurements  and  predictions  of 
blockade  radii  for  exciting  on-resonance  without  ODT.  One 
can  estimate  rb  by  equating  the  Rydberg-Rydberg  interaction 
strength  with  the  laser  excitation  bandwidth  [18,19,25,26]. 
Predictions  in  Fig.  4  are  based  on  van  der  Waals  interactions 
of  the  form  V[nt  =  C6/|R|6,  where  R  is  the  Rydberg-Rydberg 
separation  vector.  The  coefficient  Ce  is  dependent  on  the 
interacting  Rydberg  states,  | as  well  as  (for  D 
states)  the  angle  0  between  R  and  the  quantization  axis,  z.  Our 
calculations  use  interaction  strengths  from  Ref.  [19],  where  for 
D  states  we  use  C&  for  our  experimental  excitation  conditions: 
0  =  tc / 2  and,  for  z  light  polarizations,  mj  =  ±1/2.  Power¬ 
broadening  dominates  the  excitation  bandwidths,  which  are 
therefore  taken  to  be  each  state’s  two-photon  Rabi  frequency 
(ranging  from  500  kHz  for  745*  to  1.5  MHz  for  57 D). 
Agreement  between  these  predictions  and  our  measurements 
is  quite  good. 

We  further  measure  the  blockade  radius  when  exciting 
from  the  ODT.  Results  for  595*1/2  and  57 D 5/2  are  shown  in 
Fig.  5  as  a  function  of  A,  the  laser  detuning  from  the  trap-free 
transition.  As  the  detuning  A  increases,  Rydberg  atoms  are 
excited  deeper  in  the  ODT.  The  data  show  an  initial  drop  in 
rb  of  several  micrometers  with  the  application  of  the  ODT. 
With  increasing  A,  rb  decreases  for  595i/2  and  increases  for 
57D5/2.  The  trends  of  595i/2  and  51Ds/2  also  differ  near  the 
maximum  light  shift  in  the  trap  (sm2iX  ~  30  MHz  for  our  ODT 
parameters;  see  Fig.  5).  For  595i/2  we  observe  that  Rydberg 
excitation  is  possible  at  detunings  exceeding  vmax  by  about  10 
MHz  [Figs.  3(e)  and  5]. 

We  model  the  system  to  predict  the  blockade  radius, 
accounting  for  spatially  dependent  light  shifts  of  the  Rydberg 
transition.  We  consider  the  dimension  with  strongly  varying 
light  shifts  [a;  see  Fig.  1(a)].  Our  model  uses  ODT  parameters 
as  described  for  the  experiment  and  has  no  free  parameters. 


FIG.  5.  (Color  online)  Blockade  radius  rb  as  a  function  of  A, 
the  laser  detuning  from  the  trap-free  transition.  Measurements  are 
for  exciting  states  595i/2  [run  1  (±),  run  2  (x)]  and  51D5/2  (*). 
Measurements  at  A  =  0  are  without  ODT  as  in  Fig.  4;  other 
measurements  are  with  ODT.  Uncertainties  are  ±10%.  Predictions 
are  lines,  using  the  same  excitation  bandwidths  as  in  Fig.  4.  The  light 
shift  of  the  ground-Rydberg  transition  at  the  center  of  the  ODT  is 
Smax  ^  30  MHz. 


The  light  shifts  of  the  two  Rydberg  atoms,  s\  and  S2,  are 
determined  by  the  dynamic  polarizabilities  of  both  the  55i/2 
ground  state  and  the  Rydberg  state  for  1064-nm  laser  light. 
For  the  ground  state  this  is  [27]  a^s  =  47r  60  x  71 1  Aa q  ;  for  the 
Rydberg  state  it  is  calculated  via  the  ponderomotive  effect  to  be 
[28]  aRyd  =  — ■ 47T€o  x  545&q.  The  resonant  atomic  excitation 
frequencies  are  then  fi  =  fo  +  s\  and  likewise  /2  =  /o  ±  S2 
[with  /o  defined  in  Fig.  2(d)].  The  interaction  energy  between 
Rydberg  atoms  is  modeled  as 


Vint 


c6 

C Xl  -  x2)6  ’ 


(l) 


where  x\  and  X2  are  the  positions  of  the  Rydberg  atoms 
measured  transverse  to  the  ODT  beam  axis  [see  Fig.  6(a)]. 
For  5  states  Ce  >  0,  while  for  D  states  Ce  <  0  [19,24]. 

We  first  consider  5-state  atoms  for  the  case  A  <  smax.  We 
assume  a  stepwise  excitation  of  two  Rydberg  atoms,  with  the 
first  Rydberg  atom  resonantly  excited  (s  1  =  A)  at  a  position 
xi  in  the  ODT.  To  excite  a  second  atom  at  location  X2,  energy 
conservation  requires 


A  =  ^2  +  VhA/h±S,  (2) 


as  illustrated  in  Fig.  6(b),  with  an  uncertainty  8  given  by  the 
excitation  bandwidth.  The  position  X2  of  the  second  Rydberg 
atom  is  thus  determined  by  balancing  the  positive  interaction 
energy  Vint  with  a  smaller  (positive)  light  shift,  S2  <  s  1,  to 
yield  a  second  on-resonant  transition.  Among  the  solutions  for 
|xi  —  X2I,  the  smallest  one  will  manifest  itself  as  a  blockade 
radius,  so  we  take  this  one  as  the  predicted  rb.  The  results  of 
the  model,  shown  as  lines  in  Fig.  5,  agree  quite  well  with  the 
595i/2  data.  The  model  accords  with  intuition:  The  smallest 
blockade  radii  are  expected  near  the  inflection  point  of  the 
potential,  where  the  light  shift  varies  most  rapidly  with  x. 

As  noted,  for  5  states  we  measure  that  excitation  is 
possible  at  much  larger  detunings,  with  A  >  vmax.  We  model 
this  domain  by  direct  two-photon  excitation  of  two  Rydberg 
atoms;  the  off-resonant  intermediate  state  carries  one  Rydberg 
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FIG.  6.  (a)  Spatially  dependent  light  shifts  of  the  ground  and 
Rydberg  states,  with  x,  s i,  and  s2,  defined  in  text.  Panels  (b)-(d): 
Level  diagrams  for  excitation  of  two  Rydberg  atoms  from  the  ODT 
with  the  first  step  (b)  on-resonant  (/i  =  /)  for  S  states  and  [(c),  (d)] 
off-resonant  (/i  <  /)  for  (c)  S  states  and  (d)  D  states. 


excitation  [see  Fig.  6(c)].  The  intermediate- state  detuning, 
\s\  —  A |,  is  minimized  by  placing  the  first  Rydberg  atom 
at  the  center  of  the  ODT.  Energy  conservation  for  the 
direct  two-photon  excitation  of  the  Rydberg  atom  pair  then 
requires 


*1+^2  +  VM/h±8  =  2A.  (3) 

With  A  >  vmax  and  Vint  >  0,  increasing  A  gives  solutions  for 
Eq.  (3)  with  increasing  Vint  (and  thus  diminishing  r^).  Blockade 
radii  predicted  in  this  region  agree  well  with  data  in  Fig.  5.  The 
described  two-photon  excitation  process  is  similar  to  that  seen 
earlier  for  excitation  without  ODT  [Fig.  2(d)].  The  two-photon 
excitation  leads  to  the  production  of  Rydberg  atom  pairs  at 
well-defined  separations,  giving  rise  to  ring  structures  in  the 
pair  correlations.  We  can  estimate  the  relative  importance  of 
stepwise  and  two-photon  excitation  based  on  the  visibility  of 
these  structures.  Two-photon  excitation  gradually  increases  in 
importance  as  A  approaches  vmax,  so  the  rings  in  Figs.  3(a)- 
3(d)  smoothly  shrink  in  radius  and  increase  in  visibility.  The 
cusp  in  the  solid  curve  in  Fig.  5  is  an  artifact  due  to  the  model’s 
piecewise  nature,  assuming  that  two-photon  excitation  occurs 
only  when  the  more  efficient  stepwise  excitation  is  impossible. 


For  D  states  (Vint  <  0),  the  two-photon  process  described 
by  Eq.  (3)  cannot  occur  for  A  >  vmax,  as  illustrated  in  Fig.  6(d). 
Thus  the  excitation  spectra  of  D  lines  in  the  ODT  (not  shown) 
do  not  extend  beyond  A  =  vmax,  which  is  why  Fig.  5  does 
not  show  D  data  for  A  >  vmax.  However,  we  believe  that 
two-photon  processes  are  responsible  for  the  low  r ^  values 
measured  in  the  range  5  MHz  <  A  <  12  MHz.  In  this  range, 
the  laser  is  red-detuned  from  the  Rydberg  transitions  of  most 
atoms  in  the  ODT  (i.e.,  A  <  s).  Two-photon  excitation  is 
energetically  allowed  for  pairs  of  such  atoms,  because  a 
positive  value  of  s\  +  s2  —  2  A  can  be  balanced  with  a  negative 
Vint.  This  situation  can  be  described  by  a  diagram  similar  to 
Fig.  6(d)  but  with  /  <  /i,/2.  The  two-photon  process  causes 
the  Vb  data  to  be  smaller  than  predicted  by  the  stepwise  model. 
With  increasing  A  (range  12  MHz  <  A  <  vmax)  the  number  of 
atoms  that  are  available  to  this  process  (i.e.,  atoms  with  s  >  A) 
diminishes  and  the  likelihood  of  direct  two-photon  excitation 
decreases.  The  data  then  approach  the  stepwise  model  (dashed 
line  in  Fig.  5).  We  note  that  excitation  of  a  D- state  atom  at  the 
center  of  the  trap  leads  to  a  perfect  blockade,  as  neither  Eq.  (2) 
nor  Eq.  (3)  can  be  solved  for  x\  =  0  and  Vint  <  0. 

Our  one-dimensional  model  in  v  quantitatively  explains  our 
observations;  we  believe  it  captures  the  essential  physics.  In 
the  experiment  there  is  a  perpendicular  axis,  y,  along  which 
there  is  no  substantial  light-shift  gradient.  Along  this  direction 
one  might  expect  the  trap-free  blockade  radius,  but  no  clear 
anisotropy  is  observed  in  the  experimental  two-dimensional 
pair  correlations. 

In  summary,  we  have  measured  the  blockade  radius  for 
rubidium  S 1/2  and  D5/2  Rydberg  states  and  measured  the  effect 
of  laser  detunings  and  of  light-shift  potentials  on  the  Rydberg 
pair  correlation.  We  have  observed  a  variation  of  the  blockade 
radius  as  well  as  the  appearance  of  additional  structures  in  the 
pair-correlation  functions,  which  we  explain  to  be  a  result  of 
direct  two-photon  excitation  of  Rydberg  atom  pairs  at  certain 
separations.  Future  studies  may  reveal  the  dynamics  of  these 
structures  under  the  influence  of  decoherence  and  controlled 
shaping  of  the  excitation  pulses  [10,29]. 

This  work  was  supported  by  the  NSF  (PHY- 1205 5 5 9)  and 
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Production  and  trapping  of  cold  circular  Rydberg  atoms 
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Cold  circular  Rydberg  atoms  are  produced  and  magnetically  trapped.  The  trap  is  characterized  by  direct  spatial 
imaging  of  ion  distributions,  ion  counting,  and  state-selective  field  ionization.  At  room  temperature,  we  observe 
about  70%  of  the  trapped  atoms  remaining  after  6  ms.  We  measure  a  trap  oscillation  frequency  increase  of  the 
circular  Rydberg  atom  trap  relative  to  the  ground- state  atom  trap  due  to  the  larger  magnetic  moment  of  the 
circular  Rydberg  atoms.  Simulations  of  the  center-of-mass  and  internal-state  evolution  of  circular  states  in  our 
magnetic  trap  are  performed  and  results  are  in  good  agreement  with  experimental  observations. 

DOI:  10.1 103/PhysRevA.88.03 1401  PACS  number(s):  32.10.Ee,  32.10.Dk,  37.10.Jk 


Circular  Rydberg  states  [1]  correspond  to  highly  excited, 
Bohr-like  electron  trajectories  with  maximal  orbital  angular 
momenta.  Circular- state  (CS)  atoms  exhibit  a  unique  combina¬ 
tion  of  extraordinary  properties:  long  lifetimes,  large  magnetic 
moments,  and  no  first-order  Stark  shift.  These  features  have 
made  CS  atoms  attractive  for  cavity  quantum  electrodynamics 
experiments  [2,3],  where  the  CS  atoms  provide  a  near-perfect 
two-level  system  for  probing  atom-cavity  interactions.  Preci¬ 
sion  measurements  [4,5]  benefit  from  the  reduced  nuclear  and 
QED  perturbations  of  these  high-angular-momentum  states. 
The  large  angular  momenta  also  make  CS  atoms  intriguing 
candidates  for  studies  of  interactions  between  Rydberg  and 
ground-state  atoms,  such  as  the  generation  of  Rydberg-ground 
molecules  with  GHz  binding  energies  [6,7]. 

Several  applications  would  benefit  from  trapping  CS 
Rydberg  atoms.  While  optical,  electrostatic,  and  high-gradient 
magnetic  trapping  methods  have  been  considered  [8-10]  and 
demonstrated  [1 1-14]  for  atoms  in  other  Rydberg  states,  cold, 
pure  CS  atoms  have,  to  our  knowledge,  neither  been  prepared 
nor  trapped.  Additionally,  magnetic  Rydberg  trapping  has 
previously  only  been  achieved  for  mixtures  of  Rydberg  states, 
not  for  well-defined  quantum  states.  In  this  Rapid  Commu¬ 
nication  we  describe  the  generation  and  trapping  of  cold, 
pure  CS  atoms  in  a  conventional  magnetic  trap.  We  produce 
CS  atoms  at  temperatures  several  orders  of  magnitude  lower 
than  in  previous  experiments  (see,  for  example,  [1,3,15,16]) 
and  magnetically  trap  them.  We  employ  the  trapping  force 
to  measure  their  magnetic  moments,  which  are  the  largest 
used  in  atom  trapping  to  date.  We  characterize  the  trap 
using  ion  detection  schemes  that  provide  information  on  the 
center-of-mass  (COM)  oscillations  in  the  trap,  trap  losses,  and 
the  atomic  internal- state  evolution. 

We  begin  with  a  sample  of  ~107  87  Rb  atoms  magnetically 
trapped  in  the  \F  =  2,  mp  =  2)  ground  state  at  a  temperature 
~100  pK  and  density  ~1010  cm-3.  The  Ioffe-Pritchard-type 
magnetic  trap  is  generated  using  a  Z-shaped  wire  [17]. 
The  cigar-shaped  trap  has  trap  frequencies  (cox/,(Oy ,coy)  = 
2tt  x  (7.9  =b  0.3,39.6  =b  0.1,37.4)  Hz,  where  cox >  and  coy  are 
measured  and  ooz>  is  calculated.  The  trap  field  minimum  is 
measured  to  be  5.7  d=  0.6  G  and  is  located  3.5  d=  0.2  mm  from 
the  surface  of  the  Z  wire.  A  schematic  of  the  experiment  and 
a  plot  of  the  magnetic  field  strength  are  shown  in  Fig.  1.  The 
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magnetic  trapping  region  is  enclosed  by  an  electrode  structure 
consisting  of  six  individually  addressable  electrodes.  Electric 
fields  are  controlled  with  an  accuracy  of  ^30  mV/cm  in  all 
coordinates  using  voltage-to-field  calibration  factors  and  field 
zeros  obtained  from  Stark  spectroscopy  [18]  on  85 D  Rydberg 
states.  Electrode  “T”  is  used  for  Rydberg-atom  field  ionization 
and  electrode  “B”  produces  an  electric-field  ramp  required  for 
the  CS  production.  A  microchannel  plate  (MCP)  is  located 
10  cm  from  the  excitation  region,  outside  of  the  electrode 
package,  and  is  used  for  ion  detection  (30%  efficiency)  and 
imaging. 

We  generate  CS  atoms  with  principal  quantum  number 
n  =  57  by  the  crossed- fields  method  [19].  We  first  optically 
excite,  in  the  presence  of  an  electric  field  E  =  —  Ez  from 
the  B  electrode  (see  Fig.  1),  into  the  highest-lying  state 

of  the  manifold  of  n  =  57  states  that  exhibit  the  linear 

Stark  effect  (“hydrogenic  states”).  The  excitation  is  done 
with  counterpropagating  480  and  780  nm  laser  beams.  The 
fixed-frequency  780  nm  laser  is  1  GHz  off-resonance  from 
the  55*1/2  to  5^3/2  transition;  the  480  nm  laser  is  tunable  to 
excite  Rydberg  levels.  The  combined  excitation  bandwidth 
is  Vl  ~  2  MHz.  The  480  nm  beam  is  focused  to  a  waist 

of  wo  =  50  pm  and  the  780  nm  beam  is  collimated  to 

wo  =  1.5  mm,  with  powers  of  35  and  2  mW,  respectively. 
Both  beams  are  linearly  polarized  perpendicular  to  the  axis  of 
the  applied  electric  field.  At  a  field  of  ~  1.8  V/cm  the  60  R 
states  intersect  with  the  n  =  57  hydrogenic  states.  This  field 
also  mixes  enough  S  and  D  character  into  the  P  Rydberg  states 
to  enable  two-photon  excitation  from  the  5S\/2  ground  state. 
Experimental  and  calculated  two-photon  excitation  spectra  are 
shown  in  Fig.  2.  In  the  experimental  spectrum  the  hydrogenic 
states  appear  only  in  regions  where  they  intersect  with  60  P 
states.  Due  to  the  coupling  between  the  states,  the  intersections 
are  anticrossings  with  splittings  of  20  MHz  or  more,  as  seen 
in  Fig.  2(b).  For  our  initial  step  in  the  CS  production  we  excite 
to  regions  labeled  I  and  II  in  Fig.  2(a). 

After  a  15  ps  excitation  pulse  to  region  I,  the  electric  field 
is  ramped  down  over  40  ps  from  its  initial  value  of  1 .76  V /cm 
to  0  V/cm,  within  the  residual  electric-field  uncertainty. 
This  is  slow  enough  to  ensure  adiabatic  evolution  of  the 
Rydberg  atoms  within  the  anticrossings  in  Fig.  2.  Regions 
labeled  IV  are  not  chosen  because  of  state  crossings  at  lower 
electric-field  values  that  would  cause  state  mixing  during  the 
electric-field  ramp.  The  5.7  G  magnetic  field  at  the  bottom 
of  the  atom  trap  is  always  on  and  is  perpendicular  to  the 
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FIG.  1.  (Color  online)  Left:  schematic  of  the  experiment  (to 
scale).  Six  gold-coated,  copper  electrodes  enclose  a  magnetically 
trapped  atom  sample  to  control  the  electric  fields  at  the  excitation 
location.  An  MCP  is  used  for  ion  detection.  The  location  of  the  Z 
wire  and  the  counterpropagating  480  and  780  nm  beams  for  excitation 
to  Rydberg  states  are  also  shown.  Right:  plot  of  the  magnetic  field 
strength  in  the  xy  plane  through  the  center  of  the  trap.  The  linear  gray 
scale  ranges  from  \B\  =  5.5  G  (white)  to  9.5  G  (black)  in  steps  of 
0.5  G.  The  hatched  region  indicates  \B\  ^  10  G. 


ramped  electric  field.  At  the  end  of  the  ramp,  the  Larmor 
frequency  (col  =  27T  x  8  MHz)  is  larger  than  the  residual 
Stark  frequency  (cos  =  <  27r  x  3  MHz,  where  e  and 

uq  are  the  electron  charge  and  Bohr  radius,  respectively).  In 
the  adiabatic  sweep  from  the  Stark-dominated  (cos  >  col)  to 
the  Zeeman-dominated  regime  (cos  <  coL)  the  atoms  initially 
prepared  into  region  I  in  Fig.  2(a)  evolve  with  near-unity 
probability  into  the  \n  =  57,  i  =  56,  mi  =  56,  ms  =  +0.5) 
CS.  The  CS  has  good  quantum  numbers  mi  and  ms  because  in 
high- 1  states  the  Zeeman  effect  dominates  the  fine  structure. 
The  magnetic  field  stabilizes  the  CS  by  lifting  its  degeneracy 
with  any  other  hydrogenic  state.  The  angular  momentum  of 
the  trapped  CS  atom  is  parallel  to  the  magnetic  field  at  the 
trap  minimum,  which  points  along  x'  in  Fig.  1.  This  CS  is 
low-field  seeking;  the  complementary  \n  =  51,  l  =  56,  mi  = 
—56,  ms  =  —0.5)  CS  is  unsuitable  for  magnetic  atom  trapping 
because  it  is  high-field  seeking. 

Due  to  their  large  magnetic  moments,  CS  atoms  in  a 
magnetic  trap  experience  significantly  deeper  confining  poten¬ 
tials  than  ground-state  atoms.  The  magnetically  trapped  |F  = 
2,  mF  =2)  atoms  (Lande  factor  gF  =  have  a  magnetic 
moment  of  one  Bohr  magneton,  /zqs  =  ~^b-  The  CS  atoms 
described  above  have  a  magnetic  moment  /zcs  =  —  /z#(m£  + 
gsms)  =  —51/jLb,  where  gs  ~  2.  In  a  harmonic  magnetic 
trap,  the  CS  atom  experiences  a  potential  increase,  Ucs  = 
(— )Uqs  =  57£/gs»  and  a  corresponding  trap  frequency 

Mgs 

increase,  a>cs  =  VSIcoqs,  where  Ucs  and  Uqs  are  the  CS  and 
ground-state  magnetic  potentials,  and  coqs  and  coqs  are  the  CS 
and  ground- state  trap  frequencies,  respectively. 

The  minimum  of  the  ground-state  trap  is  substantially 
shifted  due  to  gravity,  while  in  the  much  stronger  CS  trap 
gravity  can  be  ignored.  The  displacement  vector  pointing  from 
the  minimum  of  the  ground- state  trap  to  the  minimum  of  the 
CS  trap  is  calculated  to  be  8r  =  ( 8x,8y,8z )  =  (710,290,0)  /xm 
(using  calculated  force  constants  and  the  12  ±  0.2  degree  tilt  of 


FIG.  2.  (a)  Experimental  Stark  spectrum  showing  the  intersection 
of  60P  with  the  n  =  51  hydrogenic  manifold.  Plotted  are  average 
counts  per  excitation  pulse  [linear  scale  ranging  from  0  (white)  to 
>3.5  (black)]  as  a  function  of  applied  electric  field  and  frequency 
offset  of  the  480  nm  laser  from  an  arbitrary  reference  frequency.  Each 
data  point  is  an  average  of  50  excitation  pulses.  The  electric-field 
axis  is  scaled  by  a  factor  of  0.98  to  match  the  calculation  in  (b). 
(b)  Energy  levels  (solid  lines)  and  excitation  rates  per  atom  calculated 
for  our  laser  polarizations  and  intensities  (overlay;  linear  gray  scale 
from  0  to  150  s-1).  In  the  excitation-rate  calculation  we  assume 
Gaussian  distributions  for  the  excitation  frequency  and  electric  field, 
with  FWHM  of  5  MHz  and  50  mV / cm,  respectively.  Minor  deviations 
in  frequency  between  (a)  and  (b)  are  attributed  to  long-term  excitation- 
laser  drifts. 


the  trap  with  respect  to  gravity;  see  Fig.  1).  This  displacement 
allows  us  to  initiate  sloshing-mode  oscillations  of  trapped  CS 
atoms  in  the  x'  and  y'  directions. 

To  observe  these  oscillations  we  first  generate  the 
| CS  —  1)  =  \n  =  57,  i  =  56,  mi  =  55,  ms  =  +0.5)  near  CS 
by  initial  excitation  to  region  II  in  Fig.  2.  We  excite  to  II 
because  it  provides  a  stronger  experimental  signal  than  I;  both 
regions  result  in  high-l  states  that  exhibit  quantitatively  similar 
trapping  behavior.  The  expected  trap  period  for  the  near  CS 
is  7cs-i  =  H=y,  which  equals  16.9  and  3.4  ms  in  x'  and 
y',  respectively.  Since  the  excitation  and  circularization  time 
interval  of  55  /zs  is  much  smaller  than  7cs-i,  the  atoms  are 
prepared  on  the  side  of  their  trapping  potential  Uqs-u  at  the 
location  —  8r.  A  field-ionization  pulse  is  applied  at  varying 
time  delays  td  from  the  end  of  the  |CS  —  1)  production.  The 
ion  positions  are  projected  onto  the  MCP,  where  they  generate 
spatially  resolved  light  pulses  on  a  phosphor  screen,  which 
is  imaged  by  a  CCD  camera.  The  COM  coordinates  of  the 
ion  distributions  in  the  MCP  plane  are  obtained  from  sums 
of  2000  images  at  each  td  [20].  Plots  of  the  x  and  y  COM 
coordinates  as  a  function  of  td  are  shown  in  Fig.  3.  The  y 
COM  data  [Fig.  3(a)]  is  fit  to  a  sinusoid  out  to  4  ms  from 
which  we  obtain  an  oscillation  period  of  3.4  ±  0.1  ms.  This  is 
in  excellent  agreement  with  the  predicted  value.  The  y  COM 
data  is  only  considered  out  to  4  ms  because  of  distortions 
in  the  ion  trajectories  thereafter,  as  suggested  by  the  x  COM 
data.  Since  8x  ~  2.58y,  we  expect  a  larger  COM  oscillation 
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FIG.  3.  (Color  online)  (a)  Experimental  COM  displacement  in  y 
vs  td  (circles)  in  the  MCP  plane  (left  axis)  with  a  fit  to  a  sinusoid  out 
to  4  ms  (solid  curve).  The  RMS  spread  of  the  data  about  the  fit  is 
0.05  mm,  indicated  by  the  representative  error  bar.  The  simulated 
COM  displacement  in  y  vs  td  in  the  object  plane  (right  axis)  is 
shown  out  to  10  ms  (dashed  curve).  A  linear  offset  is  subtracted  to 
match  the  experimental  data,  accounting  for  ion-imaging  abberations. 
(b)  Experimental  COM  displacement  in  x  vs  td.  (c)  Fraction  of 
detected  atoms  remaining  vs  td.  Based  on  the  variance  of  the  data 
points,  we  estimate  an  uncertainty  shown  by  the  representative  error 
bar. 

amplitude  in  x.  The  x  COM  coordinate  plotted  in  Fig.  3(b) 
develops  a  large  displacement  over  about  4  ms.  At  O  =  4  ms 
there  is  an  apparent  reversal  in  the  direction  of  motion  of  the  x 
COM  coordinate,  which  occurs  about  a  factor  of  2  earlier  than 
expected.  We  believe  that  the  x  COM  oscillation  amplitude 
is  so  large  that  it  brings  the  ions  produced  by  field  ionization 
far  enough  off  of  the  ion-imaging  axis  to  cause  distortion  and 
obstruction  [see  Fig.  3(c)  and  discussion].  This  prevents  a 
direct  measurement  of  the  x  COM  oscillation  period,  as  well 
as  the  y  COM  oscillation  period  after  4  ms. 

To  model  the  system,  we  perform  Monte  Carlo  simulations 
of  the  internal- state  evolution  of  a  sample  initially  in  the 
state  |CS  —  1)  in  a  T  =  300  K  radiation  background;  atom 
trajectories  inside  the  magnetic  trapping  field  are  computed 
classically.  The  internal- state  evolution  and  trajectories  are 
coupled.  The  simulations  account  for  all  radiative  transitions, 
including  thermal  photoionization,  but  no  collisions.  Typically, 
the  rate  of  thermally  induced  transitions  per  atom  is  6  x 
103  s-1.  (For  comparison,  in  a  T  =  4  K  cryogenic  enclosure 
it  would  be  0.1  x  103  s-1.)  The  rapid  state  evolution  from  the 
|CS  —  1)  due  to  the  300  K  thermal  radiation  results  in  n,  l,  and 
m  redistribution  of  the  initially  prepared  sample,  leading  to  a 
time  dependence  of  the  average  magnetic  dipole  force  as  well 
as  a  diffusive  component  of  that  force.  The  change  in  average 
m  amounts  to  less  than  1%  over  the  10  ms  observation  time. 
The  corresponding  change  in  average  oscillation  frequency 
is  too  small  to  be  seen  in  the  experiment.  The  diffusive 
component  of  the  m  evolution  is  more  significant,  reaching 
an  rms  spread  of  ±15%  about  the  average  m  over  10  ms. 
This  is  large  enough  that  we  would  expect  some  damping 
effect  on  the  COM  oscillations.  From  the  simulations  we 
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obtain  COM  oscillation  periods  and  decay  constants  of  16.5 
and  36.2  ms  in  x'  and  3.4  and  10.2  ms  in  y ',  respectively. 
In  Fig.  3(a)  we  plot  the  simulated  y  COM  oscillation  with  a 
linear  offset  subtracted  to  match  the  experiment,  accounting  for 
ion-imaging  abberations.  We  test  the  importance  of  thermally 
driven  m- state  diffusion  by  turning  off  all  radiative  transitions 
in  the  simulation,  leaving  anharmonicity-induced  dephasing 
as  the  sole  source  of  damping.  We  find  that  m  diffusion  and 
trap  anharmonicity  contribute  in  approximately  equal  parts  to 
the  decay  of  the  COM  oscillations  over  10  ms. 

To  measure  the  magnetic  trapping  time  we  plot  the  number 
of  detected  Rydberg  atoms  as  a  function  of  td.  Figure  3(c) 
shows  the  fraction  of  Rydberg  atoms  detected  relative  to  td  =  0. 
We  observe  as  much  as  70%  of  the  initially  prepared  sample 
remaining  in  the  trap  after  6  ms.  If  any  60 P  atoms  remained 
after  the  circularization  sequence,  our  simulations  show  that 
only  a~5  x  10-5  fraction  of  them  would  survive  after  6  ms. 
From  this  we  conclude  that  our  circularization  procedure  is 
very  efficient.  Between  7  and  13  ms  we  measure  a  rapid  drop 
in  signal,  as  well  as  a  significant  resurgence  at  about  16  ms. 
We  believe  that  the  rapid  drop  is  due  to  atoms  moving  out  of 
the  ion-imaging  field  of  view  (in  the  x  direction),  while  the 
resurgence  occurs  because  of  atoms  moving  back  into  view. 
This  scenario  is  consistent  with  both  the  expected  x  oscillation 
period  (16.5  ms)  as  well  as  the  early  reversal  of  the  x  COM  in 
Fig.  3(b). 

From  Fig.  3(c)  we  estimate  an  overall  loss  rate  of  >50  s-1 
for  the  trapped  CS  population.  This  loss  rate  cannot  be  ex¬ 
plained  at  a  quantitative  level,  as  low-velocity  CS  atoms  were 
not  available  until  now  and  data  on  their  collision  cross  sections 
are  not  available.  However,  the  loss  rate  is  consistent  with 
collisions  between  Rydberg  and  ground- state  atoms  being  the 
primary  loss  mechanism.  This  is  determined  from  collisional 
cross  sections  between  Rb  low -i  Rydberg  and  ground- state 
atoms  provided  in  [21]  and  the  estimated  background  Rb 
pressure.  Losses  from  Rydberg-Rydberg  collisions  are  not 
expected  at  the  low  Rydberg  atom  densities  used  (initially 
a  few  106  cm-3).  Our  simulations  of  the  state  evolution 
show  a  loss  rate  of  about  0.5  s-1  (at  300  K  and  accounting 
for  all  noncollisional  loss  channels),  which  is  negligible  in 
comparison  to  the  measured  loss  rate. 

Given  the  influence  of  the  internal- state  dynamics  on 
the  trapping  behavior,  we  probe  the  internal-state  evolution 
by  performing  state-selective  field  ionization  (SSFI)  on  the 
trapped  Rydberg  atom  sample.  A  near-linear  SSFI  ramp  is 
applied  reaching  a  field  of  >120  V/cm  over  the  course  of 
50  /xs.  Different  Rydberg  states  ionize  at  different  points  on  the 
ramp,  and  can  be  distinguished  by  arrival  time  at  the  detector. 
Figure  4(a)  shows  SSFI  traces  for  excitation  to  region  III  in 
Fig.  2,  where  only  6OP1/2  atoms  are  excited.  In  this  case, 
the  electric-field  ramp  in  the  circularization  procedure  fails  to 
populate  any  atoms  in  the  CS;  all  atoms  remain  in  6OP1/2. 
Figure  4(b)  shows  SSFI  traces  for  excitation  to  region  I  in 
Fig.  2,  where  the  CS  generation  procedure  is  expected  to 
circularize  most  of  the  atoms.  In  both  cases  we  obtain  SSFI 
traces  for  td  —  0  and  1  ms.  We  first  compare  the  integrals  of 
the  SSFI  curves.  Excitation  to  region  I  produces  a  long-lived 
signal  (0.59  remaining  fraction),  while  excitation  to  III  does 
not  (0.09  remaining  fraction),  reaffirming  the  high  yield  of 
long-lived  CS  atoms. 
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FIG.  4.  SSFI  traces  for  an  initial  excitation  into  two  different 
regions  in  Fig.  2:  (a)  into  III,  the  6OP1/2  state  and  (b)  into  I,  the  first 
anticrossing  at  the  60P -manifold  intersection  to  generate  the  CS  after 
the  circularization  procedure.  In  both  (a)  and  (b)  the  thin  and  thick 
traces  correspond  to  td  =  0  and  1  ms,  respectively.  Each  trace  is  a 
sum  of  20  000  experiments.  The  data  are  normalized  such  that  for 
td  =  0  ms  the  curves  integrate  to  the  same  value.  The  dashed  curve 
in  (a)  shows  the  field-ionization  ramp  used  in  both  experiments. 

Next  we  compare  the  structures  of  the  SSFI  curves.  In 
both  Figs.  4(a)  and  4(b),  with  td  =  0,  there  are  two  distinct 
peaks:  an  early  (22  /xs)  peak,  which  we  attribute  to  adiabatic 
field-ionization  events,  and  a  later  (29  /xs)  diabatic,  high- 
field-ionizing  peak  [22].  In  Fig.  4(b),  the  diabatic  peak  for 
td  =  1  ms  arrives  earlier  than  for  td  =  0ms,  and  is  substantially 
broadened.  This  indicates  a  reduced  average  ionization  field 
for  the  trapped  sample  at  later  times  (corresponding  to  a 
higher  average  n)  and  rapid  state  diffusion,  respectively.  An 
additional  peak  emerges  at  an  even  earlier  arrival  time  (18  /xs). 
This  peak  is  attributed  to  atoms  accumulating  in  extremely 
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high-n  Rydberg  states.  The  trends  in  Fig.  4  are  in  qualitative 
agreement  with  our  simulations,  which  also  show  a  thermally 
driven  redistribution  of  population  to  longer-lived,  very  high- 
n  states.  For  experimental  reasons  the  field  ionization  is 
performed  in  crossed  electric  and  magnetic  fields;  a  detailed 
modeling  of  the  SSFI  curves  in  Fig.  4  is  beyond  the  scope  of 
this  paper. 

In  summary,  we  have  demonstrated  a  room-temperature 
magnetic  trap  for  cold  CS  Rydberg  atoms  with  n  =  57.  The 
trap  oscillation  period  and  trapping  time  have  been  measured 
and  are  in  good  agreement  with  predicted  values  and  detailed 
simulations.  The  internal- state  evolution  of  the  atom  sample 
was  probed,  showing  the  effect  of  the  300  K  blackbody 
radiation  background  on  the  n  distribution  in  the  trap.  Our 
trapping  method  could  be  employed  in  a  4  K  cryogenic 
environment  to  suppress  thermal  transitions  and  provide  a 
source  of  long-lived  CS  atoms  for  precision  measurements 
[4,5].  Our  work  on  cold,  trapped  CS  Rydberg  atoms  paves 
the  way  for  studies  of  exotic  quantum  states.  For  example,  at 
sufficiently  high  ground- state  atom  density  high-l  Rydberg- 
ground  molecules  [6]  may  be  produced  in  the  generation  of 
the  CS.  Due  to  their  increased  mass,  these  molecules  would 
manifest  as  trailing  signals  in  the  COM  oscillations  of  the  CS  in 
the  magnetic  trap.  Another  intriguing  quantum  object,  which 
one  could  prepare  in  a  setup  similar  to  ours,  would  consist  of  a 
hollow  electron  wave  function  (that  of  the  CS  Rydberg  atom) 
surrounding  the  ionic  Rydberg-atom  core  and  a  Bose-Einstein 
condensate.  This  system  would  present  an  ideal  platform  to 
study  interactions  between  ions  and  condensates  [23-25]. 
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We  use  direct  spatial  imaging  of  cold  85 Rb  Rydberg  atom  clouds  to  measure  the  Rydberg-Rydberg 
correlation  function.  The  results  are  in  qualitative  agreement  with  theoretical  predictions  [F.  Robicheaux 
and  J.  V.  Hernandez,  Phys.  Rev.  A  72,  063403  (2005)].  We  determine  the  blockade  radius  for  states 
44£>5/2,  6OD5/2,  and  10D5/2  and  investigate  the  dependence  of  the  correlation  behavior  on  excitation 
conditions  and  detection  delay.  Experimental  data  hint  at  the  existence  of  long-range  order. 
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When  cold  atoms  are  laser-excited  to  Rydberg  states,  the 
strong  interactions  between  Rydberg  excitations  can  lead 
to  complex  many-body  entanglement  in  the  system.  The 
excitation  process  results  in  spatial  correlations  between 
excitations.  In  most  cases,  Rydberg  excitation  positions  are 
anticorrelated  such  that  no  two  excitations  are  within  a 
“blockade  radius”  rb  of  each  other.  Such  systems  are  well- 
suited  for  studying  basic  many-body  physics  [1-3],  as  well 
as  for  future  technological  applications  such  as  quantum 
computation  [4,5].  Theory  for  the  laser  excitation  process 
that  leads  to  the  Rydberg  blockade  is  considered  in 
Refs.  [6,7].  Previous  experiments  have  shown  the  blockade 
effect  to  cause  a  suppression  of  excitation  in  an  atom 
sample  [8]  and  sub-Poissonian  excitation-number  statistics 
[9].  The  underlying  energy  level  shifts  were  measured 
spectroscopically  [10]. 

We  are  presently  interested  in  the  spatial  aspect  of  the 
blockade.  In  spatially  resolved  studies,  the  distance  [11] 
and  angle  dependence  [12]  of  Rydberg-Rydberg  interac¬ 
tions  were  examined.  The  blockade  effect  was  shown  to 
be  effective  between  spatially  separated  atom  pairs  in 
adjacent  dipole  traps  [13,14].  In  this  Letter,  we  directly 
image  a  system  containing  multiple  Rydberg  excitations 
to  obtain  the  Rydberg-Rydberg  correlation  function.  We 
find  evidence  of  the  Rydberg  excitation  blockade  and 
measure  the  blockade  radius  for  several  principal  quantum 
numbers  n. 

The  experimental  setup  is  sketched  in  Fig.  1.  The  key 
feature  enabling  spatial  resolution  of  the  Rydberg  sample 
is  the  ionization  electrode:  a  beryllium-copper  needle  with 
a  rounded  tip  of  diameter  125  jmm.  This  tip  imaging  probe 
(TIP)  is  surrounded  by  a  closed  cage  of  electrodes  whose 
voltages  can  be  controlled  to  a  precision  of  several  milli¬ 
volts.  These  electrodes  allow  independent  control  of  the 
electric  field  in  three  orthogonal  directions.  To  provide 
temporal  and  spatial  electric-field  stability,  the  apertures 
in  the  electrodes  are  covered  with  wire  mesh,  which  has 
0.51  mm  grid  spacing  and  88%  transparency. 

The  experiment  runs  at  a  repetition  rate  of  30  Hz.  We 
collect 85 Rb  atoms  for  about  30  ms  in  a  magneto-optic  trap 
(MOT)  located  —500  jmm  above  the  TIP.  After  turning  off 


the  MOT  light,  we  optically  pump  the  ground  state  atoms 
into  the  5Si/2\F  =  3,  mF  =  3)  state,  with  the  quantization 
axis  along  the  beam  direction.  We  then  excite  Rydberg 
atoms  in  a  two-step  process,  using  counterpropagating  480 
and  780  nm  beams.  The  excitation  region  is  —360  jmm 
above  the  TIP.  The  beams  are  <x+  polarized  and  resonantly 
drive  the  transitions  5Si/2\F  =  3,  mF  =  3)  — ►  5 P3/2 
| F'  =4,m'F  =  4)^  nD5/2\F"  =  5 ,m''=  5).  The  red 
beam  has  a  Gaussian  beam  parameter  of  w  =  0.75  mm, 
while  the  blue  beam  is  focused  to  w0  ^  8  jmm  with  a 
theoretical  confocal  parameter  of  0.3  mm.  Because  of  the 
much  larger  size  of  the  red  beam,  the  intermediate- state 
population  is  uniform  across  the  width  of  the  excitation 
region.  We  establish  with  CCD  imaging  that  the  intensity 
distribution  of  the  blue  focus  is  only  weakly  altered  by 
diffraction  from  the  mesh.  Laser  linewidths  are  on  the 
order  of  1  MHz.  Typical  Rydberg  atom  densities  are 
— 109  cm-3,  dependent  on  atomic  state  and  laser  power. 


TIP 


FIG.  1  (color  online),  (a)  Overview  of  excitation  region  and 
ion  detection,  (b)  Exploded  view  of  electrodes  that  surround  the 
excitation  region,  used  for  electric-field  control.  The  purple, 
blue,  and  yellow  electrodes  allow  electric-field  control  in  the 
z,  y,  and  v  directions,  respectively.  Green  and  red  electrodes  are 
grounded. 
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After  excitation,  the  TIP  is  switched  to  high  voltage  to 
produce  a  strong,  radially  divergent  electric  field.  This  field 
ionizes  the  Rydberg  atoms,  and  the  ions  are  extracted  along 
the  divergent  field  lines.  The  ions  are  detected  approxi¬ 
mately  15  cm  away  by  a  microchannel  plate  detector 
(MCP),  which  has  a  spatial  resolution  of  ^  10  lines/mm. 
The  method  resembles  field  ion  microscopy  and  leads  to 
magnified  ion  images  of  the  initial  Rydberg  atom  locations. 
This  image  is  captured  by  a  CCD  video  camera,  synchro¬ 
nized  with  the  experiment. 

In  most  cases  we  use  100  ns  excitation  pulses  with  a 
subsequent  negligible  delay  before  the  field  ionization 
pulse.  Under  these  conditions  the  atoms  can  be  considered 
“frozen”  in  space  during  the  experiment.  These  parame¬ 
ters  are  later  varied  to  study  their  effect  on  the  Rydberg- 
Rydberg  correlation  function. 

To  observe  an  excitation  blockade,  it  is  essential  to  avoid 
electric  fields  during  the  excitation  pulse.  Electric  fields 
would  cause  undesired  modifications  of  the  Rydberg- 
Rydberg  interaction  potentials.  Furthermore,  electric-field 
inhomogeneity  would  induce  inhomogeneous  line  broad¬ 
ening  and  alter  the  many-body  energy  level  spectrum, 
potentially  breaking  the  excitation  blockade.  By  measuring 
the  Stark  spectra  of  nD  manifolds  as  a  function  of  the 
potentials  on  the  TIP  and  the  electrodes  shown  in  Fig.  1,  we 
zero  the  electric  field  to  below  —100  mV/cm.  This  is  low 
enough  to  avoid  the  aforementioned  problems  [15]. 

The  system’s  magnification  factor  is  determined  by  the 
combined  ion-lensing  effects  of  the  TIP,  pulsed  to  about 
400  V,  and  the  MCP  front  plate,  held  at  -800  V.  To 
calibrate  the  magnification,  we  scan  the  focused  480  nm 
beam  laterally  across  the  MCP  field  of  view  by  using  a 
mirror  with  a  piezo  actuator  on  the  horizontal  axis.  Using 
the  known  physical  displacement  of  the  beam,  we  deter¬ 
mine  the  scale  in  the  MCP  picture.  We  find  a  linear 
magnification  factor  of  45  X  between  the  excitation  region 
and  the  MCP  front  surface.  1  jmm  in  the  excitation  region 
corresponds  to  2.2  pixels  in  the  digital  CCD  images,  with 
an  uncertainty  of  10%. 

The  procedure  of  data  acquisition  and  processing  begins 
with  taking  10000  CCD  images  of  the  MCP  phosphor 
screen,  each  image  showing  the  ion  impact  positions  for 
a  single  experimental  cycle  [as  in  Fig.  2(a)].  A  CCD  image 
consists  of  these  blips  plus  background  noise  from  the 
camera.  To  eliminate  the  background,  we  determine  the 
maximum  background  level  and  subtract  it  from  the  image, 
setting  all  negatives  to  zero.  This  procedure  yields  filtered 
images,  denoted  by  X  [Fig.  2(b)].  The  autocorrelation  of 
each  filtered  image  is  calculated  and  is  denoted  A.  The 
individual  autocorrelations  A  are  then  added  to  construct  A 
[Fig.  2(d)].  We  can  write  this  as  follows,  by  using  sub¬ 
scripts  for  the  pixel  coordinates  and  superscript  a  to  enu¬ 
merate  the  images: 

A*]  =  /(«',  j,  N,  ,  (1) 

n,m 
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FIG.  2.  Stages  in  the  image  analysis  process,  using  state 
44D5/2.  All  have  the  same  scale,  (a)  Raw  image  (cropped), 
(b)  X,  filtered  image,  (c)  sum  of  filtered  images,  (d)  A,  sum  of 
autocorrelations,  (e)  B ,  autocorrelation  of  the  sum  of  images, 
(f)  A1,  normalized  autocorrelation. 


A 


i>j 


(2) 


where 


NM 

(N  -  i)(M  -  j )  9 


(3) 


N  X  M  is  the  image  size,  and  (i,  j )  is  the  displacement. 
The  normalization  factor  /  eliminates  finite-array  effects 
in  Eqs.  (2)  and  (4),  as  it  causes  flat  images  to  have  flat 
autocorrelation  functions  rather  than  pyramidal  ones. 

The  signal  A  has  three  main  structures:  a  central  peak 
due  to  the  correlation  of  a  blip  with  itself,  an  overall  cigar 
shape  due  to  the  geometry  of  the  excitation  region,  and 
structure  due  to  correlations  between  ion  positions.  It  is 
this  third  structure  that  we  wish  to  isolate  in  our  analysis. 
The  raw  autocorrelation  A  is  useful  for  qualitative  analysis, 
as  it  is  the  easiest  to  interpret  by  eye.  However,  for  quanti¬ 
tative  analysis  we  must  eliminate  the  shape  of  the  beam.  To 
do  so  we  first  sum  the  images  [Fig.  2(c)]  and  take  the 
autocorrelation  [Fig.  2(e)]: 

BiJ  =  fU,  j,  K  M)XCZXn,n)CZXn-i,m-j\  (4> 

The  autocorrelation  B  does  not  contain  any  information 
about  Rydberg-Rydberg  correlations  but  still  shows  the 
overall  shape  of  the  excitation  region.  This  allows  us  to 
divide  out  the  beam  shape  by  defining  a  normalized  A1: 
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A{j  =  AiJ/NpBiJ,  (5) 

where  Np  is  the  number  of  pictures  in  the  data  set.  The 
normalization  of  A1  is  such  that  a  value  >1  (<  1)  should 
indicate  a  correlation  (anticorrelation)  of  ion  positions. 

We  take  data  at  several  combinations  of  red  and  blue 
laser  powers.  For  the  intensity  of  the  lower  transition  beam 
we  use  intensities  of  1.7/sat  and  3.4/sat,  where  the  saturation 
intensity  7sat  =1.6  mW/cm2.  For  each  of  these  intensities 
we  measure  the  atom  number  as  a  function  of  the  intensity 
of  the  upper  transition  beam  and  observe  saturation  behav¬ 
ior  as  in  Ref.  [8].  To  acquire  images  for  the  autocorrelation 
analysis  described  above,  we  select  several  intensities  of 
the  upper  transition  beam  around  and  beyond  the  knee  of 
the  measured  saturation  curves.  We  find  that  the  autocor¬ 
relations  nearest  the  knee  tend  to  show  the  clearest 
correlation-induced  structures. 

Figure  3  shows  A  measured  for  states  44Z)5/2,  60 D5/2, 
and  70D5/2.  For  all  three  states,  we  observe  a  ring-shaped 
dip  around  the  central  self-correlation  peak  in  the  autocor¬ 
relation.  Such  dips  are  present  in  most  of  our  data,  with  the 
degree  of  visibility  depending  on  excitation  parameters. 
The  absence  of  dips  in  some  data  indicates  that  they  are  not 
artifacts  of  the  imaging  or  image  processing.  We  find  that 
laser  frequency  fluctuations  can  make  the  dips  change  in 
depth  and  diameter.  Furthermore,  laser  power  affects  the 
results,  as  described  below. 

To  quantitatively  compare  the  autocorrelations  A1  of 
different  states,  in  Fig.  4  we  plot  7(r),  defined  as  the 
angular  average  of  A1.  7(r)  is  analogous  to  the  radial 
autocorrelation  function  calculated  in  Ref.  [6].  The  self¬ 
term  fills  the  region  r  <  5  jmm .  We  do  not  remove  the  self¬ 
term  due  to  the  risk  of  producing  a  false  blockade  signal. 
Between  5  and  10  jmm,  each  curve  has  a  dip  indicating  an 
anticorrelation  in  Rydberg  atom  positions,  as  was  qualita¬ 
tively  observed  in  Fig.  3.  The  minima  of  7(r)  do  not 
approach  zero  in  part  because  we  are  projecting  a  three- 
dimensional  sample  onto  the  detector  plane,  which  reduces 
the  visibility  for  even  a  perfect  blockade. 

Figure  4  shows  that  the  blockade  radius  increases  with 
principal  quantum  number  n.  However,  there  is  no  obvious 
choice  for  a  method  to  measure  a  blockade  radius  from 

44  D5/2  60D5/2  70D5/2 


■  10|Lim 

FIG.  3.  Autocorrelations  A  for  447)5/2,  60D5/2,  and  70D5/2 .  In 
the  data  sets  shown  we  use  a  laser  intensity  of  3.47sat  for  the 
lower  transition  and  the  lowest  intensity  setting  (nearest  the 
saturation  knee)  for  the  upper  transition. 


each  curve,  as  the  curve’s  shape  is  determined  by  both  the 
blockade-induced  dip  and  the  self-term. 

To  determine  a  good  measure  for  the  blockade  radius, 
we  simulate  the  excitation-blockade  process  for  our  ge¬ 
ometry  and  z-projection  conditions,  by  using  a  stochastic, 
nonquantum  model.  We  find  that  the  inflection  point  after 
the  minimum  of  each  7(r)  curve  is  close  to  the  input  hard- 
sphere  blockade  radius.  We  thus  use  this  inflection  point  as 
our  measure  of  the  blockade  radius.  We  note  that  in  the 
simulations  the  self-term  diminishes  the  dip  depth  and  can 
shift  the  dip  position.  If  the  input  blockade  radius  is  close 
to  the  self-term  width,  the  radius  indicated  by  the  inflection 
point  is  1-2  jmm  higher  than  the  blockade  radius.  When  the 
blockade  radius  is  much  larger  than  the  self-term  width,  the 
inflection  point  method  underestimates  the  input  blockade 
radius  by  ^  1  jmm. 

For  the  states  44D5/2,  60 D5/2,  and  70D5/2  we  measured 
the  radii  on  nine,  six,  and  nine  curves.  One  additional  curve 
was  discarded  for  70D5/2;  although  its  dip  minimum  was 
consistent  with  other  data,  the  position  of  its  inflection 
point  was  unclear.  Five  additional  curves  for  44D5/2 
were  thrown  out  because  there  was  no  dip,  as  discussed 
below.  The  measured  blockade  radii  are  shown  in  Fig.  5. 
The  error  bars  indicate  the  10%  uncertainty  from  the 
magnification  calibration;  this  systematic  uncertainty 
dominates  the  uncertainty  due  to  measurement  statistics 
by  a  factor  of  at  least  2. 

In  Fig.  5,  we  compare  the  measurements  to  predictions 
based  on  an  excitation  bandwidth  S  vL  ranging  from  5  to 
12  MHz  and  interaction  strengths  as  determined  from 
Ref.  [15].  The  blockade  radius  for  the  van  der  Waals 
interaction  has  been  calculated  to  be 

rAW(2)(n*)1111/6 

r*  “  J  •  <6) 

where  n*  =  n  —  St  is  the  effective  principal  quantum 
number,  St  is  the  quantum  defect,  and  A  W®  is  a  scaled, 
second-order,  state-dependent  van  der  Waals  coefficient 


FIG.  4  (color  online).  Angular  average  of  A1,  7(r),  for  447)5/2, 
6OD5/2,  and  707)5/2.  Curves  correspond  to  images  in  Fig.  3. 
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FIG.  5  (color  online).  Comparison  of  blockade  radius  mea¬ 
surements  (  +  )  with  predictions  (X)  based  on  Ref.  [15]. 


defined  in  Ref.  [15].  A  W®  is  scaled  such  that  it  would  be 
independent  of  n*  in  the  absence  of  resonant  level  shifts. 
The  assumed  value  for  S  vL  has  little  effect  due  to  the  sixth- 
root  dependence. 

Figure  5  shows  that  the  measured  blockade  radii  are 
within  a  factor  of  2  of  the  predictions.  Further,  the  amount 
by  which  the  blockade  radius  increases  between  the  probed 
n  levels  is  as  expected.  The  measured  blockade  radii 
exceed  the  calculated  radii  by  2-3  jmm .  The  difference 
may  be  because  Eq.  (6)  models  the  interaction  between 
two  atoms  and  excludes  many-body  effects  [3].  Ates  et  al. 
have  previously  seen  evidence  that  theory  based  on  pair¬ 
wise  interactions  underestimates  the  blockade  radius  [7]. 

To  test  if  atomic  motion  or  the  degree  of  excitation 
saturation  modifies  the  observed  correlation  behavior,  we 
vary  laser  powers,  excitation  pulse  duration,  and  field 
extraction  delays.  In  general,  increasing  the  laser  power, 
or  the  excitation  duration  at  constant  power,  diminishes  the 
dip  depth  (defined  as  the  curve  asymptote  minus  the  dip 
minimum).  For  instance,  increasing  the  excitation  duration 
for  state  70D5/2  from  100  to  400  ns  causes  the  dip  depth  to 
decrease  from  0.24  to  0.04.  This  behavior  could  be  caused 
by  an  increased  z  depth  of  the  saturated  part  of  the  excita¬ 
tion  region.  In  our  simulations  the  degree  to  which  the  dip 
depth  is  diminished  depends  on  the  assumed  beam  focus 
quality.  Adding  wings  to  a  Gaussian  focus,  as  suggested  by 
Fig.  2(c),  yields  better  qualitative  agreement  with  our  data. 

Laser  power  has  a  greater  effect  for  the  case  of  excitation 
to  44Z)5/2  than  for  the  other  states:  In  five  7(r)  curves  the 
dip  disappeared  entirely.  We  expect  that  this  is  partly  due  to 
saturating  the  atomic  transition  for  44Z)5/2  more  strongly 
than  for  60 D5/2  or  70D5/2.  It  is  also  possible  that,  of  the 
three  states  we  studied,  44D5/2  is  the  most  susceptible  to 
laser  parameters  due  to  having  the  shallowest  dip,  closest 
to  the  radius  of  the  self-term.  Another  possibility  is  that  the 
blockade  efficiency  is  diminished  due  to  the  proximity  of 
the  Forster  resonance,  as  suggested  by  Refs.  [3,16]. 


Multiple  rings  Maxima  showing  angular  structure 

FIG.  6.  Some  autocorrelations  (A)  showing  possible  long- 
range  and  angular  structure  for  state  44D5/2. 


By  varying  the  delay  between  excitation  and  field  ion¬ 
ization,  we  test  the  longevity  of  the  correlations.  They  are 
remarkably  long-lived.  The  dip  depth  decreases  by  about 
half  over  10  juls,  despite  possible  atomic  motion  and  ion¬ 
ization  during  this  time  [17,18].  This  result  awaits  a  future 
theoretical  explanation. 

Some  of  our  data,  using  linear  laser  polarization  and  a 
few  megahertz  of  detuning,  show  additional  structures.  The 
autocorrelation  in  Fig.  6(a)  exhibits  multiple  rings,  indicat¬ 
ing  long-range  Rydberg-Rydberg  correlations.  The  auto¬ 
correlation  in  Fig.  6(b)  hints  at  multiple  maxima  in  a  ring 
about  the  center,  indicating  possible  angular  dependence  of 
the  Rydberg-Rydberg  correlation  function.  There  have 
been  predictions  of  crystal-like  structure  in  Rydberg 
systems  [19],  but  the  detailed  excitation  conditions  are 
important  in  generating  these  structures.  At  present  we 
are  working  to  improve  our  experimental  control  in  order 
to  systematically  study  these  effects. 

In  summary,  we  have  obtained  spatially  resolved  images 
of  blockaded  Rydberg  atom  samples  and  evaluated  the 
Rydberg-Rydberg  correlation  function.  Having  observed 
a  blockade,  we  measured  the  blockade  radius  for  several 
quantum  states.  We  have  also  tested  the  dependence  of  the 
blockade  on  parameters  such  as  laser  power,  excitation 
duration,  and  detection  delay.  Future  work  could  include 
investigation  of  longer-range  correlation  structures, 
extending  beyond  one  blockade  radius,  as  in  Fig.  6.  One 
could  further  attempt  to  excite  a  Rydberg  crystal  as  pro¬ 
posed  in  Ref.  [19]. 
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We  study  angular-momentum  couplings  in  87Rb2  Rydberg  molecules  formed  between  a  Rydberg  and  55i/2 
ground-state  atom.  We  use  a  Fermi  model  that  includes  5-wave  and  P-wave  singlet  and  triplet  scattering  of 
the  Rydberg  electron  with  the  55i/2  atom,  along  with  the  fine- structure  coupling  of  the  Rydberg  atom  and 
hyperhne- structure  coupling  of  the  55i/2  atom.  We  discuss  the  effects  of  these  couplings  on  the  adiabatic 
molecular  potentials.  We  obtain  bound-state  energies,  lifetimes,  and  electric  and  magnetic  dipole  moments  for 
the  vibrational  ground  states  of  the  87Rb (nD  +  55i/2)  molecules  in  all  adiabatic  potentials,  with  fine  and  hyperfme 
structure  included.  The  hyperhne  structure  gives  rise  to  mixed  singlet- triplet  potentials.  We  also  study  the  effect 
of  the  hyperhne  structure  on  the  deep  35-wave-  and  3P -wave-dominated  adiabatic  molecular  potentials,  which 
support  high-f  87Rb2  Rydberg  molecules. 
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I.  INTRODUCTION 

Rydberg  molecules  formed  by  the  low-energy  scattering 
of  a  Rydberg  electron  and  a  ground-state  atom  constitute  a 
distinct  class  of  molecular  states  that  have  in  recent  years 
become  the  subject  of  significant  theoretical  and  experimental 
interest.  Developments  include  studies  of  diatomic  Rydberg 
molecules  in  low- angular-momentum  Rydberg  S  states  [1,2], 
P  states  [3],  and  D  states  [4,5],  the  realization  of  coherent 
bonding  and  dissociation  of  5-type  molecules  [6],  and  the 
first  observation  of  a  permanent  electric  dipole  moment  in 
a  homonuclear  molecule  [7].  Polyatomic  Rydberg  molecules 
have  also  been  generated  [8]  and  employed  in  a  demonstration 
of  the  continuous  transition  between  a  few-body  to  many-body 
regime  in  an  ultracold  quantum  gas  [9] . 

The  theoretical  framework  for  these  molecules  is  generally 
well  established.  The  interaction  between  a  low-energy  Ryd¬ 
berg  electron  and  ground-state  atom  can  be  described  using  a 
Fermi  pseudopotential  approach  [10-12].  In  the  Fermi  model, 
the  ground-state  atom  is  treated  as  a  8 -function  perturber  of 
the  Rydberg-electron  wave  function,  resulting  in  oscillatory 
potential  curves  with  localized  minima  capable  of  sustaining 
bound  molecular  states.  In  alkali  systems,  the  potential  curves 
are  strongly  affected  by  low-energy  35-wave  and  3P-wave 
electron-atom  scattering  resonances  [13-15].  The  influence 
of  these  scattering  resonances  on  the  potentials  has  been 
studied  in  Rb2.  The  3S  interaction  leads  to  so-called  trilobite 
molecules  [12],  which  are  very  long  range  (their  size  is  on  the 
order  of  n2,  where  n  is  the  principal  quantum  number  of  the 
Rydberg  level).  The  3P  interaction  produces  potentials  [16] 
that  are  about  an  order  of  magnitude  deeper  and  about 
a  factor  of  five  shorter  range  than  the  trilobite  potentials. 
The  molecular  potentials,  wave  functions,  and  electric  dipole 
moments  for  both  Rb2  and  Cs2  have  also  been  calculated  using 
a  Green’s  function  approach,  accounting  for  the  finite  size  of 
the  perturbing  ground-state  atom  via  an  effective  short-range 
electron-atom  interaction  potential  [17,18]. 

In  addition  to  the  electron-atom  scattering  interaction,  the 
molecular  potentials  and  properties  of  long-range  Rydberg 
molecules  are  dependent  on  the  Rydberg-atom  wave  function 
as  well  as  the  angular-momentum  couplings  of  the  Rydberg- 
and  ground- state  constituents.  Rydberg  molecules  exhibit 


a  broad  range  of  different  angular  momentum  coupling 
schemes.  For  low -i  Rydberg  molecules  (t  <  2  in  rubidium), 
the  angular-momentum  coupling  configurations  span  three 
Hund’s  cases  [(a)-(c)],  dictated  by  the  relative  strength  of  the 
Rydberg  atom’s  fine- structure  coupling  compared  to  that  of 
the  e~  +  55i/2  scattering  interaction.  The  Rb (nDj  +  5S\/2) 
molecules  are  unique  among  the  low -l  molecules  because 
their  Rydberg  fine- structure  couplings  are  comparable  to 
the  scattering  interaction  strength,  and  they  fall  anywhere 
between  two  Hund’s  cases  (a)  and  (c)  by  a  mere  change  in 
principal  quantum  number  n.  For  high-£  Rydberg  molecules, 
the  molecular  binding  interaction  is  stronger  than  the  fine- 
structure  coupling,  and  is  comparable  to  the  hyperhne  coupling 
of  the  ground- state  perturber.  Inclusion  of  the  5S\/2  hyperhne 
coupling  in  the  model  generates  additional  adiabatic  potentials 
of  mixed  triplet  and  singlet  character.  For  both  high-  and 
low -l  molecules,  the  hyperhne  structure  results  in  additional 
adiabatic  potentials  deep  enough  to  sustain  bound  states. 

In  this  work  we  present  a  study  on  the  inhuence  of  angular- 
momentum  couplings  on  the  properties  of  long-range  Rydberg 
molecules.  We  hrst  describe  a  Fermi  model  for  Rydberg 
molecules  with  the  relevant  angular-momentum  couplings 
included.  In  our  analysis  we  include  singlet  and  triplet  5- 
and  P-wave  scattering.  By  selectively  enabling  the  different 
interactions  in  the  model,  the  effects  of  each  of  these  individual 
interactions  on  the  adiabatic  molecular  potentials  are  revealed. 
Typical  vibrational- state  wave  functions,  binding  energies, 
lifetimes,  and  dipole  moments  are  discussed  with  an  emphasis 
on  the  role  of  angular-momentum  couplings.  We  conclude  with 
describing  hyperhne  effects  on  the  relatively  deeply  bound 
trilobite  Rydberg  molecules  [12]. 


II.  FERMI  MODEL 

We  describe  diatomic  Rydberg  molecules  with  a  Fermi 
model  [  1 0, 1 1  ]  taking  into  account  the  angular-momentum  cou¬ 
plings  in  the  Rydberg  atom  and  perturbing  ground-state  atom 
system  whose  strengths  are  comparable  to  the  Rydberg  e~-\- 
perturber  interaction.  A  schematic  of  the  relevant  couplings  is 
shown  in  Fig.  1(a).  The  perturbing  87Rb  5S\/2  atom  is  located 
at  a  position  Z  from  the  ionic  core  of  the  Rydberg  atom, 
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FIG.  1.  (a)  Angular  momentum  coupling  scheme  for  diatomic 
87Rb (nDj  +  5Si/2,F)  Rydberg  molecules.  The  relevant  interactions 
are  circled.  Here,  AO  and  A1  denote  e~  +  55i/2  scattering  interactions 
involving  the  mn  =  0  (5-wave  and  P-wave)  and  the  \ntn  \  =  1  (P- 
wave  only)  components  of  the  Rydberg  electron’s  state,  respectively, 
and  H  denotes  the  hyperfine  interaction  of  the  55i/2  atom.  FS 
denotes  the  fine- structure  coupling,  (b)  States  in  the  mk  =  mj\  + 
ms 2  +  mi2  =  +1/2  subspace  and  their  relevant  interactions.  In  the 
left  column,  horizontal  gray  bars  are  placed  between  (m7im52)  or 
(m52mj2)  for  states  in  neighboring  rows  that  are  coupled  by  either  the 
scattering  or  the  hyperfine  interaction.  The  Xs  in  the  right  column 
indicate  the  interactions  that  have  terms  diagonal  in  (rrijim^m^). 


which  is  fixed  at  the  origin.  The  internuclear  axis  is  along  z. 
The  orbital  and  spin  angular  momenta  of  the  Rydberg  atom  are 
denoted  by  L  i  and  S\ ,  respectively.  For  L  i  <  2  and  within  the  n 
range  of  interest,  the  Rydberg-atom  fine  structure  is  of  the  same 
order  as  the  e~  +  5S\/2  scattering  interaction  and  is  therefore 
included.  The  hyperfine  coupling  of  the  electron  spin  S2  and 
nuclear  spin  I2  of  the  perturbing  5S\/2  ground-state  atom  is 
also  included,  because  it  is  several  GHz  and  is,  in  most  cases, 
stronger  than  the  Rydberg  electron’s  fine- structure  coupling 
and  the  e~  +  5S\/2  scattering  interaction.  The  orbital  angular 
momentum  of  the  5S\/2  atom  is  L2  =  0.  The  Rydberg  atom’s 
hyperfine  structure  decreases  as  n~ 3;  for  the  lowest  5  states 
relevant  to  our  work  it  does  not  exceed  several  MHz  [19,20], 
and  it  is  much  lower  for  higher-^  states  [21].  The  Rydberg 
atom’s  hyperfine  structure  is  therefore  not  included. 

For  a  Rydberg  atom  with  its  ionic  core  at  the  origin  and  the 
Rydberg  electron  located  at  r,  and  with  a  5S\/2  atom  located 
at  Zz,  the  Hamiltonian  is  written  as 


energy-dependent  scattering  phase  shift.  For  the  calculations  in 
the  present  work  we  use  nonrelativistic  scattering  phase  shifts 
8i= 0  and  <$/=i  provided  by  I.  I.  Fabrikant  based  on  Ref.  [17]. 
The  electron  momentum  is  given  by  k  =  y/—\/n q2  +  2/r 
(atomic  units)  in  the  classically  allowed  range  of  the  Rydberg 
electron  and  k  =  0  elsewhere.  Here,  is  the  effective  Rydberg 
quantum  number  of  the  level  of  interest.  To  account  for  config¬ 
uration  interactions,  we  employ  basis  sets  {\n,L\,J\,nij\)  0 
that  include  all  Rydberg  levels  with  effective 
quantum  numbers  \n*  —  n^\  <  2.5.  For  the  Rydberg  atom,  all 
L\,  /1,  and  mj\  are  included  (as  in  Ref.  [24]),  and  for  the 
perturber  atom  all  ms2  and  m^  are  included.  The  singlet  (5) 
and  triplet  (7)  channels  of  the  e~  +  5  S\/2  scattering  interaction 
include  projectors  I(s,t),  defined  as  1 7  =  Si  •  §2  +  f ,  which 
has  an  eigenvalue  of  one  (zero)  for  the  triplet  (singlet)  states, 
and  Is  =  1  —  It-  The  operators  S\  and  S2  are  the  spins  of 
the  Rydberg  electron  and  5S\/2  atom,  respectively,  and  I  is 
the  identity  operator.  The  last  term  in  Eq.  (1)  accounts  for  the 
hyperfine  interaction  of  the  perturber.  The  87  Rb  5  S1/2  perturber 
atom  has  nuclear  spin  I2  with  I2  =  3/2,  and  hyperfine  levels 
F<  =  1  and  F>  =  2,  with  a  hyperfine  coupling  parameter 
A  =  h  x  6.8  GHz /F>  =hx3A  GHz  (in  SI  units). 

Only  levels  with  mn  =  0,±1  components  have  nonva¬ 
nishing  wave  functions  or  wave  function  gradients  on  the 
internuclear  axis.  The  5- wave  interactions  couple  Rydberg 
states  with  mn  =  0  components.  The  P-wave  interactions 
couple  states  with  mn  =  0  components  through  the  radial 
derivative  of  the  Rydberg  wave  function,  and  states  with  mn  = 
=L  1  components  through  the  polar-angle  derivative  of  the  wave 
function.  The  electron  scattering  term  in  Eq.  (1)  conserves 
mji  +  ms 2,  while  the  hyperfine  term  conserves  ms2  +  mi2.  The 
full  Hamiltonian  conserves  :=  mj\  +  ms2  +  m* 2.  Hence, 
the  Hilbert  space  can  be  broken  up  into  subspaces  of  fixed 
quantum  number  As  an  example,  the  subspace  for 
mk  =  +1/2  and  its  couplings  via  the  5-wave,  P-wave,  and 
hyperfine  interactions  is  shown  in  Fig.  1(b).  The  gray  bars 
in  Fig.  1(b)  illustrate  that  the  couplings  are  organized  in  an 
interlaced  block-diagonal  structure  in  the  magnetic  quantum 
numbers.  Since  the  couplings  via  the  scattering  terms  require 
\mji\  ^  3/2,  the  87Rb2  Rydberg  molecules  have  |m^|  <  7/2. 


III.  ADIABATIC  POTENTIALS 


H(r,Z)  =  H0  +  2nA‘s(k)83(r  -  Zz)!-, 

i=S,T 

+  6?rA^(/fc)S3(r  -  Zz)S/  ■  VI,- 

i—S,T 

+  AS2-I2.  (1) 

In  the  unperturbed  Rydberg  Hamiltonian  Hq  we  use 
published  quantum  defects  [22],  which  account  for  core 
penetration  and  fine  structure.  For  0  5  we  use  the  fine- 
structure  correction  of  hydrogen  as  well  as  a  small  quantum 
defect  to  account  for  core  polarization,  81  =  0J5aD/l5  [22], 
with  a  dipolar  polarizability  for  Rb+  of  ao  =  9.023  atomic 
units  [23].  The  energy-dependent  5- wave  (/  =  0)  and  F-wave 
(/  =  1)  scattering  lengths  (As  and  Ap,  respectively)  have  the 
general  form  A/(k)  =  —  tan<$//k2/+1,  where  5/  is  the  /-  and 


Adiabatic  molecular  potentials  V/(Z)  (/  is  an  arbitrary  label) 
are  obtained  by  diagonalizing  the  Hamiltonian  in  Eq.  (1)  for  a 
grid  of  perturber  atom  positions,  Zz,  through  the  extent  of  the 
Rydberg  wave  function.  To  highlight  the  effects  of  the  different 
terms  in  Eq.  (1)  on  V/(Z),  in  Figs.  2  and  3  we  show  adiabatic 
potentials  calculated  for  D-type  Rydberg  molecules  with  the 
different  interaction  terms  in  Eq.  (1)  selectively  turned  on. 

First,  we  consider  adiabatic  potentials  without  hyperfine 
coupling.  Figure  2(a)  shows  adiabatic  potentials  for  the 
31D  +  55i/2  molecule  with  only  the  35  interaction  turned 
on  and  no  fine-structure  coupling.  This  results  in  three 
degenerate  oscillatory  potentials,  one  corresponding  to  each 
triplet  state,  and  a  flat  potential  corresponding  to  the  singlet 
state  (which  has  no  35  interaction).  The  triplet  potential  curves 
are  similar  to  those  calculated  in  Ref.  [12],  in  which  the 
5-wave  scattering  length  —  tan  8s/k  is  taken  to  first  order  in 
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FIG.  2.  Adiabatic  potentials  for  the  3 ID  +  55i/2  molecule  with 
the  following  interaction  terms  in  Eq.  (1)  selectively  turned  on 
(without  the  hyperfine  interaction):  (a)  3S  scattering,  (b)  3S,  lS,  3P,  and 
lP  scattering,  and  (c)  3S ,  lS,  3P,  and  lP  scattering  with  fine-structure 
coupling. 

the  electron  momentum  [11].  With  an  appropriate  choice  of 
the  zero-energy  5-wave  scattering  length,  the  3S  interaction 
reproduces  measured  binding  energies  of  vibrational  ground 
states  of  5-type  Rydberg  molecules  [1].  Similarly,  the  3S 
interaction  with  the  addition  of  the  fine  structure  reproduces 
vibrational  ground  states  of  D-type  Rydberg  molecules  [4]. 
The  effect  of  the  fine- structure  coupling  on  the  molecules  is 
discussed  further  below. 


FIG.  3.  Binding  adiabatic  potentials  for  the  87Rb(30D3/2  + 
5Si/2,  F  =  1,2)  molecules  with  fine  and  hyperfine  structure  included, 
with  3S  scattering  only  (top  row)  and  with  3S,  ]S,  3P,  and  lP  scattering 
(bottom  row).  The  hyperfine  coupling  leads  to  the  shallow  adiabatic 
potentials  of  mixed  singlet-triplet  character.  The  shallow  potentials 
are  different  for  the  F  =  1  and  F  =  2  hyperfine  levels.  The  deep 
potentials  have  pure  triplet  character  and  do  not  depend  on  the 
hyperfine  structure. 


Figure  2(b)  shows  adiabatic  potentials  resulting  from  the 
3S ,  lS ,  3P ,  and  lP  scattering  interactions  turned  on  and  no 
fine-structure  coupling.  At  smaller  internuclear  separations, 
the  influence  of  R-wave  scattering  becomes  more  significant 
due  to  higher  electron  energies  closer  to  the  Rydberg  atom’s 
ionic  core  (and  therewith  larger  wave  function  gradients).  The 
effect  of  the  3P  scattering  interaction  is  evident  in  Fig.  2(b), 
where  the  inner  wells  become  notably  deeper,  while  the 
outermost  well  remains  relatively  unaffected.  The  increasing 
contribution  of  R-wave  scattering  at  smaller  Z  generates 
deep  molecular  potentials.  The  repulsive  lS  and  attractive  lP 
scattering  interactions  turn  the  flat  singlet  potential  in  Fig.  2(a) 
into  an  oscillatory  singlet  potential  with  maxima  above  and 
wells  below  the  dissociation  threshold,  as  seen  in  Fig.  2(b). 

Figure  2(c)  shows  adiabatic  potentials  resulting  from  the 
the  addition  of  the  fine- structure  coupling  of  the  Rydberg 
atom  to  the  3S ,  lS ,  3R,  and  lP  scattering  interactions.  The 
top  and  bottom  plots  show  31  Dj  +  55*1/2  potentials  for  the 
j  =  5/2  and  3/2  fine-structure  states,  respectively.  Qualita¬ 
tively,  the  fine- structure  coupling  splits  the  molecular  bonding 
strength  (adiabatic-potential  depth)  of  the  fine-structure-free 
case  between  the  two  fine-structure  levels,  resulting  in  less 
deep  potentials.  The  splitting  ratio  depends  on  which  Hund’s 
case  is  more  relevant  (see  Sec.  V). 

The  hyperfine  interaction  of  the  ground-state  perturber  in 
Eq.  (1)  mixes  the  singlet  and  some  triplet  scattering  channels, 
resulting  in  the  replacement  of  the  pure  singlet  potentials  in 
Fig.  2  with  shallow  adiabatic  potentials  of  mixed  singlet-triplet 
character.  As  an  example,  in  Fig.  3  we  show  the  adiabatic  po¬ 
tentials  for  the  3OD3/2  +  55*1/2, R  =  1,2  molecules  including 
the  hyperfine  interaction.  The  deep  adiabatic  potentials  are 
of  pure  triplet  character  and  are  unaffected  by  the  hyperfine 
interaction  of  the  55i/2  atom.  The  F  =  1  shallow,  mixed 
singlet-triplet  potentials  are  always  deeper  than  the  F  =  2 
shallow  potentials.  The  shallow  potentials  can  typically  sustain 
a  few  bound  states  that  should  be  observable  in  experiments. 

IV.  QUASIBOUND  MOLECULAR  STATES  AND  LIFETIMES 

Due  to  the  lack  of  an  inner  potential  barrier,  the  molecular 
vibrational  states  have  the  character  of  metastable  scat¬ 
tering  resonances  in  potentials  that  are  unbound  on  the 
inside.  Figure  4  shows  the  (deep)  adiabatic  potential  for 
the  87Rb(31D3/2  +  55i/2,R  =  2)  molecule  and  its  quasibound 
states.  These  are  qualitatively  similar  to  those  of  5-type 
molecules  [8].  The  molecular  wave  functions  consist  of  low- 
amplitude  standing  waves  formed  by  outgoing  and  reflected  in¬ 
going  waves  in  the  region  Z  <  lOOO^o,  and  a  high- amplitude, 
resonantly  enhanced  portion  in  the  outer  adiabatic-potential 
wells  at  Z  ~  1500<2o.  The  latter  are  identified  with  quasi¬ 
bound  molecular  vibrational  states  that  are  metastable  against 
tunneling-induced  decay  (decay  into  the  region  Z  <  lOOCLzo). 

The  resonances  are  found  by  computing  the  wave  function 
phase  at  a  fixed  location  in  the  unbound  region  (we  use  Z  = 
300tfo)  as  a  function  of  energy  W.  The  phase  and  its  derivative 
are  plotted  as  a  function  of  energy  in  the  right  and  middle 
panels  of  Fig.  4,  respectively.  The  quasibound  molecular  states 
occur  at  energies  at  which  the  phase  undergoes  sudden  changes 
of  AO  =  7r,  and  are  centered  at  energies  where  the  derivative 
of  the  phase  is  maximal  (circles  in  the  middle  panel  of  Fig.  4). 
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FIG.  4.  (Left)  Adiabatic  potential  of  the  87Rb(31D3/2  + 
551/2, F  =  2)  Rydberg  molecule  and  vibrational  wave  functions. 
Each  wave  function  corresponds  to  a  narrow  scattering  resonance, 
characterized  by  a  sudden  change  in  the  wave  function  phase  by  tt  in 
the  unbound,  inner  region  of  the  potential.  (Right)  Wave  function 
phase  at  location  Z  =  300 a0.  (Middle)  The  maxima  of  d<&/dW, 
indicated  by  circles,  are  used  to  determine  resonance  widths  and 
lifetimes.  Several  broad  resonances  (hatched  regions)  are  spread  out 
over  the  displayed  energy  range. 


The  resonances  obey  a  single-level  Breit-Wigner  formula  with 
frequency  linewidths  of  the  quasibound  molecular  states  given 
by  ry  =  2/ (h  x  d^/dW)  [25],  corresponding  to  lifetimes  of 
t  =  (h/2)d<&/dW .  The  lifetimes  scale  as  the  Wigner  tunneling 
time  delay  [26].  Resonances  with  larger  slopes  dQ/dW  in 
Fig.  4  correspond  to  longer-lived  quasibound  states.  In  addition 
to  the  sharp  resonances  there  are  broad  resonances,  as  indicated 
by  hatched  regions  in  the  middle  panel  of  Fig.  4.  While  the 
broad  resonances  are  not  likely  to  cause  observable  effects  in 
molecular  spectra,  they  add  to  the  total  phase  change  over  an 
energy  range.  For  instance,  within  the  energy  range  displayed 
in  Fig.  4  the  broad  resonances  account  for  a  phase  change 
of  37T  and  the  quasibound  molecular  states  for  a  change  of 
llyr,  corresponding  to  a  total  change  of  1471.  The  distinction 
between  broad  and  narrow  resonances  may,  in  practice,  depend 
on  experimental  parameters  such  as  excitation  bandwidth  and 
atom  temperature  (which  affects  Frank-Condon  factors). 

In  Table  I  we  give  the  binding  energies,  linewidths,  decay 
times,  and  average  intemuclear  separation  of  the  scattering 
resonances  shown  in  Fig.  4.  Each  quasibound  state  is  assigned 
a  vibrational  quantum  number  v,  sequentially  increasing  from 
0  for  the  most  strongly  bound  state  in  the  outermost  well  to 
10  for  the  most  weakly  bound  state.  Generally,  one  expects 
bound  states  furthest  from  the  dissociation  threshold  to  have 
the  longest  lifetimes.  Here,  the  ground  and  first  excited 
states  of  the  outermost  potential  well  (labeled  v  =  0  and 
2,  respectively)  are  well  confined  and  have  lifetimes  in  the 
range  of  hours  and  seconds,  respectively.  These  lifetimes  only 
reflect  tunneling-induced  decay.  The  actual  lifetimes  of  these 
molecular  states  are,  in  fact,  much  shorter  due  to  additional 
decay  mechanisms,  such  as  radiative  decay  of  the  Rydberg 
state  and  collisions  with  ground-state  atoms  [27].  An  early 
dissociation  process  via  the  energy  exchange  between  the 


TABLE  I.  For  the  87Rb(31D3/2  +  5 Sl/2,F  =  2)  Rydberg 
molecule,  we  show  the  vibrational  quantum  number  v,  binding 
energy,  linewidth  Tv,  decay  time,  and  bond  length. 


V 

Energy 

(MHz) 

Linewidth 

(MHz) 

Decay  time 
Os) 

(Z) 

(ao) 

0 

-  100.52 

1.79X10-11 

8.89x  109 

1647.5 

1 

-83.59 

2.18x  10-2 

7.30 

1160.0 

2 

-  68.03 

6.00x  10“8 

2.65  xlO6 

1655.3 

3 

-  56.27 

2.37  xl0“2 

6.72 

1324.0 

4 

-  45.24 

1.48x10-’ 

1.08 

1413.0 

5 

-41.58 

1.87x10“’ 

0.85 

1527.9 

6 

-36.91 

1.01x10° 

0.16 

1169.3 

7 

-26.58 

1.28x10“’ 

1.24 

1428.4 

8 

-  18.48 

9.99x  10“2 

1.59 

1467.0 

9 

-9.14 

7.78x  10“2 

2.05 

1481.8 

10 

-1.48 

3.71  xl0“2 

4.29 

1697.0 

Rydberg  electron  and  ground-state  atom  has  also  been  used 
to  explain  shorter  molecular  lifetimes  observed  in  experi¬ 
ments  [28] .  Nevertheless,  the  relatively  long  lifetimes  and  large 
Frank-Condon  factors  associated  with  these  outermost  states 
make  them  the  easiest  to  isolate  experimentally  [1,4,9].  There 
are  also  several  resonances  in  the  inner  potential  wells  (v  =  1 
and  3)  with  lifetimes  on  the  order  of  that  of  the  Rydberg  atom 
(which  for 87 Rb  3ID3/2  is  about  20  /xs  in  a  300  K  black-body 
radiation  field).  The  resonances  at  higher  energies  (v  =  4  —  6) 
exhibit  shorter  lifetimes  because  of  the  smaller  potential  barrier 
through  which  they  more  readily  tunnel  inward.  Surprisingly, 
above  these  short-lived  resonances  additional  resonances 
with  longer  lifetimes  appear  (v  =  7  —  10).  Here,  the  inner 
oscillatory  wells  act  like  an  aperiodic  Bragg  reflector  of  the 
molecular  wave  functions,  resulting  in  unexpectedly  long- 
lived  resonances  near  the  dissociation  threshold.  The  lifetime 
of  these  states  is  largest  when  a  Bragg  reflection  condition 
is  met.  This  occurs  in  a  range  of  Rydberg  principal  quantum 
numbers  n  at  which  the  periodicity  of  the  vibrational  wave 
function  approximately  equals  that  of  the  Rydberg-electron 
wave  function.  This  Bragg-reflection  has  previously  been 
described  as  an  internal  quantum  reflection  process  in  5-type 
molecules  [8]. 


V.  HUND’S  CASES  FOR  nD  RYDBERG  MOLECULES 

Hund’s  coupling  cases  are  widely  used  for  classification 
of  angular  momentum  couplings  in  diatomic  molecules  [29] . 
how-l  diatomic  Rydberg  molecules  (£  <  2  in  rubidium) 
exhibit  a  variety  of  coupling  cases,  determined  by  the  relative 
strength  of  the  Rydberg  atom’s  fine-structure  coupling  to  the 
e~-\-  perturber  scattering  interaction.  For  87 Rb (nDj  +  5 £1/2) 
molecules  the  fine-structure  coupling  is  comparable  to  the 
scattering  interaction  strength.  Due  to  this,  D-type  molecules 
trend  from  Hund’s  case  (c)  at  large  n ,  where  the  fine-structure 
coupling  exceeds  the  scattering  interaction  strength,  to  Hund’s 
case  (a)  for  n  <  35,  where  the  scattering  interaction  strength 
exceeds  the  fine- structure  coupling.  In  previous  work,  we 
observed  87 Rb(nDj  +  55i/2,F  =  2)(v  =  0)  molecules  in  tran¬ 
sition  between  the  two  Hund’s  coupling  cases  (a)  and  (c)  [4].  In 
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Z  (units  of  a0)  Z  (units  of  a0) 


FIG.  5.  (Color  online)  (a)  Adiabatic  potentials  for  87Rb(22D7  + 
55i/2)  (left)  and  87Rb(40D7-  +  5Sl/2)  (right)  with  inter¬ 

actions  and  no  HFS  interaction,  (b)  Binding  energies  for  the  v  =  0 
ground  vibrational  state  of  the  nD5/2  +  5S\/2  (blue  squares)  and 
nD3/ 2  +  5Si/2  (red  triangles)  molecular  potentials  versus  n.  The  D 
fine-structure  splitting  is  also  plotted  (black  circles). 


this  section,  we  focus  on  the  potentials  and  quasibound  states 
of  the  87Rb (nDj  +  5*Si/2),  j  =  3/2  and  5/2  molecules  in  the 
two  limiting  Hund’s  cases  and  in  the  transition  regime. 

Figure  5(a)  shows  VJ-(Z)  for  the  j  =  3/2  and  5/2  (22 Dj  ± 
55'i/2)  molecules  (left)  and  j  =  3/2  and  5/2  (40  Dj  ±  55i/2) 
molecules  (right)  calculated  with  all  interaction  terms  in 
Eq.  (1),  excluding  hyperfine- structure  coupling.  At  high  n, 
the  molecules  trend  towards  Hund’s  case  (c),  where  the 
dominant  adiabatic  molecular  potentials  are  reduced  by  the 
product  of  two  Clebsch-Gordan  coefficients  of  the  type 
(mi  =  0 ,ms  =  d=l/2|  j,mj  =  ±1/2).  This  leads  to  adiabatic 
potentials  whose  depths  carry  spin-dependent  factors  1/(21  ± 
1)  for  j  =  i  —  1/2  and  (t  +  l)/(2 i  +  1)  for  j=l  +  1/2  [25]. 
For  nD  (i  =  2)  molecules  at  high  n,  the  depth  ratio  of  the 
potentials  for  j  =  3/2  and  j  =  5/2  approaches  2/3.  This  is 
seen  in  the  depths  of  the  outermost  wells  of  the  high -n  j  =  3/2 
and  5/2  (40 Dj  ±  5Si/2)  molecular  potentials  in  Fig.  5(a).  For 
decreasing  n,  the  fine-structure  splitting  increases  as  n~ 3  while 
the  scattering  interaction  strength  increases  as  n~ 6,  and  the 
molecules  tend  towards  Hund’s  case  (a).  The  j  =  3/2  and  5/2 
(22 Dj  +  5Si/2)  molecular  potentials  in  Fig.  5  exhibit  this  case, 
where  the  scattering  interaction  strength  is  large  compared  to 
the  fine-structure  splitting.  In  the  low -n  limit,  the  j  =  3/2 
potential  becomes  notably  deeper  than  the  j  =  5/2  potential. 

The  binding  energies  of  the  vibrational  ground  states  in  the 
outermost  potential  wells,  W/>u= o,  closely  track  the  depth  of 


those  wells  and  therefore  mirror  the  transition  of  the  molecule 
between  the  two  Hund’s  cases.  In  Fig.  5(b),  we  plot  W/>v= o  for 
the  j  =  3/2  and  5/2  (nDj  +  5Si/2)  molecules  for  a  range  of 
n  from  Hund’s  case  (a)  (low-n)  to  Hund’s  case  (c)  (high-n). 
In  the  high-n  limit,  the  molecular  binding  energies  W/>v= o 
for  both  fine-structure  levels  approximately  scale  as  n-6, 
inversely  with  the  atomic  volume.  In  the  low-n  limit,  the  lowest 
binding  energy  for  the  lower  fine- structure  level  (j  =  3/2) 
is  larger  and  continues  to  scale  as  ~n-6,  while  that  of  the 
upper  fine-structure  level  (j  =  5/2)  trends  towards  a  ~n~3 
scaling,  approaching  the  scaling  of  the  fine-structure  splitting. 
For  the  Rb  (nDj  +  5Si/2)  molecules,  the  binding  energies  of 
the  vibrational  ground  states  in  the  outermost  potential  wells 
for  j  =  3/2  and  j  =5/2  are  approximately  equal  at  n  =  34 
[see  Fig.  5(b)].  Further,  due  to  the  described  Hund’s  case 
behavior,  at  sufficiently  low  n  the  inner  wells  of  adiabatic 
potentials  of  the  nD^/2  ±  5Si/2  molecule  become  deep  enough 
to  support  vibrational  states  that  are  more  deeply  bound  and 
long  lived  than  the  ground  states  of  the  outermost  well  of  the 
nDs/2  +  5S\/2  potentials. 


VI.  ELECTRIC  AND  MAGNETIC  DIPOLE  MOMENTS 


Rydberg  molecules  present  the  only  known  case  of 
homonuclear  molecules  with  permanent  electric  dipole  mo¬ 
ments  [7,12].  For  high-£  Rydberg  molecules  dipole  moments 
on  the  order  of  \03eao  are  predicted  to  exist  [12].  Smaller 
permanent  dipole  moments  arise  in  low -l  S- ,  P-,  and  D- 
type  molecules  from  fractional  admixing  of  high-l  state 
character.  Permanent  electric  dipole  moments  of  ~0.5eao 
have  previously  been  measured  in  rubidium  S-type  Rydberg 
molecules  [7].  Dipolar  cesium  Rydberg  molecules  with 
electric  dipole  moments  of  ~5  —  50eao  [2]  have  also  been 
prepared.  In  this  section  we  calculate  both  the  electric  and 
magnetic  dipole  moments  for  n  D  Rydberg  molecules,  with  all 
terms  in  the  Hamiltonian  in  Eq.  (1)  included. 

We  obtain  the  adiabatic  electric  [<7;,z(Z)]  and  magnetic 
[l 'JLitZ(Z )]  dipole  moments  in  the  diagonalization  of  Eq.  (1). 
The  dipole  moments  of  a  molecular  state  v  follow  from  the 
expectation  values  of  di  z(Z)  and  /x/  z(Z)  over  the  vibrational 
wave  function  densities, 


di,v  =  j  l^i,v(Z)|2rfi,z(Z)dZ 
Rv  =  j  l^,v(Z)|2Rz(Z)JZ. 


(2) 


Electric  dipole  moments  for  the  ground  vibrational  states  of 
the  j  =  3/2  and  5/2  87Rb (nDj  +  5Si/2, F  =  1,2)  deep  (pure 
triplet)  and  shallow  (mixed  singlet- triplet)  molecular  potentials 
are  shown  in  Fig.  6.  In  Figs.  6(a)  and  6(b)  we  show  the  binding 
energies  and  representative  potentials  with  wave  functions, 
respectively,  for  all  angular-momentum  coupling  cases  that 
arise  from  Eq.  (1).  We  note  that  for  87Rb (nDj  +  5S\/2,F  = 
1,2)  molecules  the  hyperfine  quantum  numbers  F  are  well 
defined  because  the  hyperfine  coupling  is  much  larger  than 
the  molecular  binding.  As  seen  in  Figs.  6(a)  and  6(b),  the 
delineation  between  Hund’s  cases  (a)  and  (c)  is  shifted  to 
lower  n  for  the  shallow  potentials  (which  are  due  to  hyperfine- 
induced  mixing  of  singlet  and  triplet  states  at  the  5Si/2 
atom).  Essentially,  the  generally  weaker  scattering  interaction 
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FIG.  6.  (Color  online)  (a)  Binding  energies  for  v  =  0  in  the  outermost  potential  wells  vs  n,  with  all  terms  in  Eq.  (1)  included.  Symbol 
legend:  Solid:  deep,  triplet  potentials.  These  are  identical  for  F  =  1  and  F  =  2.  Hollow:  shallow,  mixed  singlet-triplet  potentials.  These  are 
different  for  F  =  1  and  F  =  2.  Black  squares:  j  =  5/2.  Red  circles:  j  =  3/2.  Large  size:  \mk\  =  1/2.  Medium  size:  \mk\  =  3/2.  Small  size: 
\mk\  =  5/2.  The  nDj  fine-structure  splitting  is  plotted  for  reference.  Numbers  indicate  degeneracies  summed  over  all  mk.  (b)  Representative 
adiabatic  potentials  and  wave  functions  for  87Rb(29D7  +  5S\/2,F  =  1,2 )(v  =  0),  for  j  =  5/2  (left  column)  and  j  =  3/2  (right  column).  The 
pure  triplet  potentials  are  the  same  for  F  =  1  and  F  =  2  (top  row),  while  the  mixed  singlet-triplet  potentials  are  generally  shallow  and  different 
for  F  =  1  (middle  row)  and  F  =  2  (bottom  row).  The  gray  bars  on  the  right  indicate  binding  energy,  to  visualize  that  the  deep  potentials  are 
closer  to  Hund’s  case  (a)  (j  =  2/2  potential  deeper  than  j  =  5/2  potential)  than  the  shallow  ones  (j  =  5/2  potentials  deeper  than  j  =  2/2 
potentials),  (c)  Electric  dipole  moments  diyV  for  v  =  0  in  the  outermost  potential  wells  vs  n ,  with  all  terms  in  Eq.  (1)  included.  Symbol  legend 
as  in  (a).  The  blue  line  through  the  data  for  the  deep  j  =  3/2  potentials  is  an  allometric  fit  with  exponent  —2.4.  (d)  Magnetic  dipole  moments 
/xijV  for  v  =  0  in  the  outermost  potential  wells  vs  n,  with  all  terms  in  Eq.  (1)  included.  Symbol  legend  as  in  (a).  We  only  show  data  for  positive 
mk  (results  for  negative  mk  are  the  same  with  flipped  sign).  There  are  no  degeneracies  in  /x/  v=0- 


associated  with  the  mixed  singlet-triplet  cases  pushes  those 
molecules  more  towards  Hund’s  case  (c). 

The  electric  dipole  moments  of  the  deep  molecular  poten¬ 
tials  are  eightfold  degenerate  when  summed  over  all  and 
those  of  the  shallow  potentials  are  fourfold  degenerate  [see 
Fig.  6(c)].  It  is  noted  that  with  decreasing  n  the  degeneracies 
in  energy  and  electric  dipole  moment  become  increasingly 


lifted.  This  may  be  attributed  to  a  stronger  configuration 
mixing  at  low  n  caused  by  the  relative  increase  of  the  e~-\- 
perturber  scattering  term  in  Eq.  (1).  The  only  case  in  which  the 
dipole  moments  of  the  87Rb(/t  Dj  +  5  S\/2 ,  F  =  1,2)  molecules 
exhibit  a  clear  scaling  behavior  in  n  is  for  the  j  =  3/2  deep 
potentials  (which  are  the  same  for  F  =  1  and  F  =  2).  A  fit 
to  the  topmost  data  set  in  Fig.  6(c)  gives  an  n~2A  scaling, 
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similar  to  a  result  found  previously  for  87Rb2  5 -type  Rydberg 
molecules  [7].  Here,  the  scattering-induced  mixing  between 
the  atomic  Rydberg  levels  gives  rise  to  a  significant  change  in 
the  electric  dipole  moments  as  a  function  of  n .  In  the  range  n  > 
30,  the  electric  dipole  moments  for  both  types  of  shallow  po¬ 
tentials  of  87Rb(ftZ)3/2  +  5Si/2,  F  =  1,2)  also  approximately 
scale  as  n~2A.  The  electric  dipole  moments  for  the  upper  fine- 
structure  component  j  =  5/2  have  less  clear  scaling  trends 
[squares  in  Fig.  6(c)].  In  particular,  the  electric  dipole  moments 
for  the  deep  87Rb(nDs/2  +  5S\/2,F  =  1,2)  potentials  do  not 
exhibit  a  clear  scaling  behavior.  This  is  likely  a  result  of  the 
transition  from  Hund’s  case  (a)  to  Hund’s  case  (c). 

The  magnetic  moments  for  the  molecular  states  are  shown 
in  Fig.  6(d).  The  magnetic  moments  are  nondegenerate  due 
to  the  different  g  factors  of  the  involved  types  of  spins.  The  n 
dependence  again  reflects  the  transition  in  angular-momentum 
coupling  behavior  between  the  two  Hund’s  cases  (a)  and  (c). 
We  expect  that  experiments  in  weak  electric  and  magnetic 
fields  can  reveal  the  electric  and  magnetic  dipole  moments  of 
the  vibrational  states.  The  dipole  moments  computed  in  this 
work  are  for  the  weak-field  limit,  i.e.,  results  are  expected  to 
be  accurate  as  long  as  molecular  Zeeman  and  Stark  shifts  are 
smaller  than  other  relevant  energy  scales  (such  as  the  energy 
splitting  between  adjacent  vibrational  states). 


VII.  HYPERFINE-STRUCTURE  EFFECTS  IN 
DEEP  35-  AND  3P-DOMINATED  POTENTIALS 

The  trilobite  molecules  [12]  generated  by  the  3S  interaction 
are  long  range  and  on  the  order  of  10  GHz  deep,  which  is 
comparable  to  the  hyperfine  interaction  of  the  55*1/2  perturber. 
It  is  therefore  of  interest  to  explore  hyperfine  effects  on 
the  trilobite  adiabatic  potentials.  In  Figs.  7(a)  and  7(b)  we 
show  these  potentials  without  and  with  hyperfine  coupling  for 
n  =  30,  for  all  relevant  values  of  rriji  +  ms2  and  m k  =  rtij\  + 
tyis2  +  m,2,  respectively.  While  the  hyperfine  interaction  does 
not  affect  the  general  shape  and  depth  of  the  trilobite  potentials, 
with  hyperfine  structure  included  there  are  three  instead  of 
only  one.  Since  the  number  of  3P  potential  curves  also  triples, 
the  crossing  pattern  between  3S  and  3P  potentials  becomes 
considerably  more  complex,  as  seen  in  Fig.  7.  Considering 
that  the  modulations  near  the  bottom  of  the  trilobite  potentials 
are  several  100  MHz  deep,  which  is  sufficient  to  support 
individual  bound  vibrational  states,  the  hyperfine  structure  is 
expected  to  have  profound  effects  on  the  detailed  vibrational 
level  structure  of  these  molecules.  In  Fig.  7(b),  one  may  expect 
to  find  long-lived  ^-dominated  states  in  the  gray  regions. 
Further,  as  indicated  by  the  circles  in  Fig.  7(b),  the  crossing 
locations  and  the  detailed  coupling  behavior  of  i  =  3  Rydberg 
levels  (which  are  optically  accessible  from  low-lying  atomic 
states)  to  molecular  states  in  the  long-range  potentials  also 
strongly  depend  on  the  hyperfine  structure  of  the  system  under 
investigation. 

In  the  present  work  we  have  not  included  the  dependence 
of  the  3P  £-+Rb(5Si/2)  scattering  phase  shift  on  J  (where 
J  =  0, 1,2)  [17].  We  have  performed  a  preliminary  calculation 
that  indicates  that  the  outermost  several  wells  of  the  adiabatic 
potentials,  which  support  the  metastable  molecular  states,  are 
not  substantially  altered  by  using  /-dependent  3P  scattering 


Z  (units  of  a0) 


FIG.  7.  (Color  online)  High-f  adiabatic  potentials  near  n  =  30  of 
87Rb2  (Rydberg  +  55i/2)  molecules  (a)  without  and  (b)  with  hyperfine 
structure  included.  The  plots  indicate  the  trilobite  potentials  [12],  the 
dominant  types  of  scattering  interactions  leading  to  deep  potentials, 
the  asmyptotic  states  of  the  potentials,  the  regions  where  bound 
trilobite  molecules  may  be  found  (gray  areas),  and  the  intersections 
between  trilobite  potentials  and  F  (l  =  3)  lines  (dashed  circles). 


phase  shifts.  We  also  find  that  the  crossings  between  the  3P 
potentials  and  the  3S  trilobite  potentials  shift  in  position  by 
amounts  up  to  about  50  ao. 


VIII.  CONCLUSION 

In  summary,  we  have  systematically  explored  the  role 
of  fine- structure  and  hyperfine- structure  couplings  in  87Rb2 
molecules  formed  between  a  Rydberg  and  55*1/2  ground- 
state  atom.  As  has  been  done  extensively  in  previous  work, 
we  have  treated  the  e~  +  5  S\/2  scattering  with  a  Fermi 
model  that  includes  5-wave  and  P-wave  singlet  and  triplet 
scattering.  The  fine  structure  mostly  influences  the  behavior 
of  low-f  Rydberg  molecules.  We  have  explored  in  detail 
how  87Rb(ftD  +  5S\/2)  molecules  realize  Hund’s  cases  (a) 
and  (c).  The  hyperfine  structure  originates  in  the  perturber 
atom  and  therefore  has  consequences  for  all  types  (low-£  and 
high-f!)  of  Rydberg  molecules.  In  the  case  of 87  Rb(nD  +  5S\/2) 
molecules,  mixing  of  singlet  and  triplet  potentials  results  in  a 
set  of  shallow  adiabatic  potentials,  whose  quasibound  states 
should  be  experimentally  observable  (in  addition  to  those  in 
the  hyperfine-independent  pure  triplet  potentials).  We  have 
obtained  electric  and  magnetic  dipole  moments,  which  could, 
in  future  work,  be  measured  spectroscopically  in  weak  electric 
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and  magnetic  fields.  The  hyperfine  structure  has  also  been  seen 
to  alter  the  deep,  ^-scattering-induced  trilobite  potentials  as 
well  as  their  crossing  behavior  with  i  =  3  Rydberg  states; 
experimental  spectroscopic  studies  should  reveal  these  details. 
Since  molecular  level  energies  and  properties  are  very  sensitive 
to  the  scattering  phase  shifts  used  in  the  Fermi  model,  we 
expect  that  spectroscopy  of  Rydberg  molecules  can  serve 
as  a  sensitive  tool  to  provide  measurement-based  input  to 
future  theoretical  studies  in  low-energy  electron  scattering. 


The  effects  of  the  J  dependence  of  the  3P  scattering  phase 
shifts  could  also  be  addressed  in  future  work. 
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We  observe  long-range  homonuclear  diatomic  nD  Rydberg  molecules  photoassociated  out  of  an 
ultracold  gas  of 87 Rb  atoms  for  34  <  n  <  40.  The  measured  ground-state  binding  energies  of  87Rb {nD  + 
55^2)  molecular  states  are  larger  than  those  of  their  87Rb {nS  +  5Si/2)  counterparts,  which  shows  the 
dependence  of  the  molecular  bond  on  the  angular  momentum  of  the  Rydberg  atom.  We  exhibit 
the  transition  of  87Rb(*rD  +  55'1/2)  molecules  from  a  molecular-binding-dominant  regime  at  low  n  to  a 
fine- structure-dominant  regime  at  high  n  [akin  to  Hund’s  cases  (a)  and  (c),  respectively].  In  the  analysis,  the 
fine  structure  of  the  nD  Rydberg  atom  and  the  hyperfine  structure  of  the  55^2  atom  are  included. 

DOI:  10.1 103/PhysRevLett.l  12.163201  PACS  numbers:  34.50.Cx,  31.10,+z,  33.80.Rv,  34.20.Cf 


Cold  atomic  systems  have  opened  new  frontiers  at  the 
interface  of  atomic  and  molecular  physics.  Of  particular 
interest  are  a  recently  discovered  class  of  long-range, 
homonuclear  Rydberg  molecules  [1,2].  Formed  via  an 
attractive  interaction  between  a  Rydberg  electron  and  a 
ground-state  atom  [1],  these  molecules  are  among  the 
largest  ever  observed  with  intemuclear  separations  of 
several  thousand  Bohr  radii.  Their  distinctive  binding 
mechanism,  which  is  unlike  conventional  covalent,  ionic, 
and  van  der  Waals  bonds  between  ground-state  atoms, 
results  in  loosely  bound  molecules  whose  properties  mimic 
those  of  their  constituent  Rydberg  atoms.  The  discovery 
of  these  molecular  bonds  has  been  likened  to  a  new 
ultracold  chemistry  [3]  and  has  spurred  a  significant 
amount  of  theoretical  [4,5]  and  experimental  interest 
[6-9].  Nondegenerate,  low  angular  momentum  Rydberg 
states  (orbital  angular  momentum  €  <  2  in  rubidium) 
produce  molecules  with  ~  MHz  binding  energies  and 
permanent  electric  dipole  moments  of  a  few  Debye.  The 
€  —  0  molecules  were  first  observed  by  photoassociation 

[10]  of  cold  87Rb  atoms  [2].  Recently,  molecular  states  with 
€  —  1  and  high  electron  energies  have  also  been  excited  via 
bound-bound  transitions  in  85Rb2  [9].  For  high- Rydberg 
states,  so-called  trilobite  molecules  with  giant  permanent 
electric  dipole  moments  of  several  kilo-Debye  and  ~  GHz 
binding  energies  are  predicted  to  exist  [1].  These  have  yet 
to  be  observed,  though  dipolar  molecules  have  been 
prepared  in  Rb  [7]  and  Cs  [8]. 

The  relevant  interaction  was  first  described  by  Fermi 

[11]  to  help  explain  pressure-induced  energy  shifts  of 
Rydberg  absorption  lines  in  a  gas  [12].  The  de  Broglie 
wavelength  of  the  Rydberg  electron  (position  r)  is  much 
larger  than  that  of  a  heavy  ground-state  atom  (position  R) 
that  lies  within  the  Rydberg  atom’s  volume,  and  their 
interaction  can  be  approximated  as  a  low-energy  v-wave 
scattering  process  (scattering  length  as).  The  interaction 
is  described  with  a  Fermi-type  pseudopotential  [1,13] 
^pseudo (r)  =  2nas S3  (r-R),  where  p- wave  and  higher-order 


scattering  [4]  are  neglected.  For  negative  as ,  the  interaction 
can  lead  to  bound  molecular  states  [1,13]. 

The  bond  strengths  and  electric  and  magnetic  moments 
of  Rydberg  molecules  are  strongly  dependent  on  the 
angular-momentum  couplings  in  the  Rydberg  and 
ground-state  atom  constituents.  In  the  present  work,  we 
focus  on  long-range  87Rb2  molecules  formed  by  an  nD 
Rydberg  and  a  5 S\/2  ground-state  atom.  Within  the  low-^f 
class  of  these  molecules,  the  nD  ones  have  the  highest 
binding  energies,  which  generally  increase  with  *?  due  to 
the  \]2£  +  1  scaling  of  the  angular  wave  functions 
YJl=0(6  =  0).  The  angular-momentum  coupling  spans 
three  Hund’s  cases  as  €  varies  from  0  to  2.  The  nS\/2  + 
5 S\/2  molecules  are  akin  to  Hund’s  case  (b),  because  they 
have  L  =  0  and  total  electron  spin  S  =  1.  The  nPj  +  5 S\/2 
molecules  are  akin  to  Hund’s  case  (c),  because  the  fine 
structure  coupling  is  larger  than  the  scattering  interaction 
energy.  The  nDj  +  5Si/2  molecules  present  a  unique 
system,  because  they  transition  from  Hund’s  case  (c)  at 
large  n  to  Hund’s  case  (a)  at  n  <  30,  where  the  scattering 
interaction  energy  becomes  increasingly  greater  than  the 
fine  structure  coupling.  The  F  and  higher-*?  states  of  Rb 
intersect  with  a  sole,  deeply  bound  (trilobite)  potential 
which  accumulates  most  of  the  level  shifts  [1],  negating 
the  existence  of  potential  wells  similar  to  those  below  the 
S ,  P,  and  D  asymptotes.  In  our  model,  we  include  the 
fine  structure  of  the  Rydberg  atom  and  the  hyperfine 
structure  of  the  5S\/2  perturber  atom,  both  of  which  have 
observable  effects. 

To  excite  nD  Rydberg  molecules,  we  first  prepare  a 
sample  of  ~105  magnetically  trapped  87  Rb  atoms  in  their 
\F  =  2,  mF  =  2)  ground  state  at  a  temperature  <  17  pK  and 
peak  density  >5  x  1011  cm-3.  Optical  excitation  to  nD 
Rydberg  states  is  accomplished  via  a  two-photon  transition 
from  the  5 Si/2  ground  state  by  using  780-  and  480-nm  laser 
beams.  The  780-nm  laser  frequency  is  fixed  ~1  GHz  off 
resonance  from  the  55,1/2  to  5 P3/2  transition,  and  the  480-nm 
laser  frequency  is  scanned  to  excite  either  Rydberg  atoms  or 
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molecules.  The  combined  excitation  bandwidth  is  #2  MHz. 
The  780-nm  laser  has  a  power  of  ~500  fi W  and  is  collimated 
to  a  full  width  at  half  maximum  (FWHM)  of  3.5  mm.  The 
480-nm  beam  has  a  power  of  ~35  mW  and  is  focused  to  a 
FWHM  of  89  ±  5  /nm  into  the  cigar-shaped  atom  sample, 
which  has  a  FWHM  diameter  of  28  /um  and  an  aspect  ratio 
of  #1:3. 

The  atom  sample  is  enclosed  by  six  individually 
addressable  electrodes  used  to  control  the  electric  field 
at  the  excitation  location,  described  in  previous  work 
[14].  We  zero  the  electric  field  to  within  <200  mV / cm  by 
Stark  spectroscopy  on  59 D  Rydberg  states  [15].  This 
ensures  that  quadratic  Stark  shifts  of  the  nD  Rydberg 
levels  in  the  n  range  of  interest  are  <2  MHz.  In  a  single 
experiment,  the  atom  sample  is  illuminated  by  2  —  3  //s 
long  laser  pulses  followed  by  electric-field  ionization 
of  Rydberg  atoms  and  molecules  [16].  The  signal  ions 
are  extracted  and  detected  by  a  microchannel  plate 
located  10  cm  away.  We  use  one  ground-state  atom 
sample  for  a  series  of  55  individual  experiments  at  a 
single  480-nm  frequency  step;  density  loss  during  one 
series  is  negligible. 

The  photoassociation  of  a  Rydberg  atom  and  ground- 
state  atom  pair  into  a  bound  molecular  state  occurs  when 
the  excitation  laser  is  detuned  from  the  atomic  Rydberg  line 
by  an  amount  equal  to  the  molecular  binding  energy.  In 
Fig.  1,  we  show  an  experimental  spectrum  in  the  vicinity 
of  the  atomic  35D5/2  Rydberg  line.  A  prominent  satellite 
line  emerges  at  —  38  =b  3  MHz,  which  is  assigned  to  the 
(: 35D5/2  ±  5Si/2)(v  =  0)  molecule,  where  v  —  0  denotes 
the  vibrational  ground  state.  This  binding  energy  is  #1.6 
times  larger  than  that  of  the  (355'1/2  +  55,1/2)(^  =  0) 
molecular  state  measured  in  previous  experiments  [2,6]. 


Frequency  offset  (MHz) 

FIG.  1.  Spectrum  centered  on  the  35D5/2  atomic  Rydberg  line 
showing  the  87Rb(35D5/2  +  5Si/2){v  =  0)  molecular  line  at 
-38  d=  3  MHz.  The  vertical  error  bars  are  the  standard  error  of 
three  sets  of  55  individual  experiments  at  each  frequency  step. 
The  error  in  the  binding  energy  is  equal  to  the  average  long-term 
frequency  drift  observed  over  one  full  scan. 


The  larger  binding  energy  reflects  the  deepening  of  the 
molecular  potential  with  f. 

A  series  of  nD5/2  Rydberg  spectra  for  34  <  n  <  42  is 
shown  in  Fig.  2(a)  (right).  The  lowermost,  redshifted  lines 
for  n  —  40  and  below  are  assigned  to  the  molecular  states. 
Molecular  lines  are  not  discernable  in  the  n  =  42  and  41 
spectra,  because  the  line  broadening  due  to  residual  fields 
and  the  laser  linewidth  exceeds  the  molecular  binding 
energies  for  these  states.  Additional  satellite  lines  are 
expected  but  are  likely  obscured  in  Figs.  1  and  2(a)  by 
the  broadening  of  the  atomic  lines  as  well  as  artificial 
signals  at  ±20  MHz  due  to  weak,  symmetric  side  peaks  in 
the  480-nm  laser  spectrum  (caused  by  a  Pound-Drever-Hall 
stabilization  loop).  Features  near  -20  MHz  are  assigned  to 
molecular  lines  only  if  they  are  significantly  stronger  than 
the  artificial  signal  at  ±20  MHz. 


Principal  quantum  number  (n) 

FIG.  2.  (a)  Right:  Spectra  centered  on  nD5/2  atomic  Rydberg 

lines  for  the  indicated  values  of  n  and  identified  molecular  lines 
(squares).  Left:  Selected  spectra  from  the  plot  on  the  right  for 
states  with  identified  molecular  lines.  Error  bars  are  obtained  as  in 
Fig.  1.  (b)  Binding  energies  obtained  from  Gaussian  fits  to  the 
molecular  lines  identified  in  (a)  vs  n.  An  allometric  fit  (solid 
curve)  to  the  experimental  binding  energies  yields  a  n~5-9±0A 
scaling.  Also  shown  are  theoretical  binding  energies  for  the 
87 Rb (nD5/2  ±  5S\/2){v  =  0)  molecular  states  with  as0  —  -I4a0 
(hollow  circles  and  dotted  curve),  (c)  Peak  number  of  detected 
ions  on  the  nD5/2  atomic  Rydberg  line  (solid  diamonds,  left  axis) 
and  ratio  of  molecular  and  atomic  line  strengths  (hollow 
diamonds,  right  axis)  vs  n. 
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Figure  2(b)  shows  the  molecular  binding  energies 
measured  in  Fig.  2(a)  vs  n.  One  may  expect  the  binding 
energies  to  be  proportional  to  the  probability  density  of 
the  Rydberg  electron  wave  function,  which  scales  as 
~n~6.  An  allometric  fit  to  the  data  in  Fig.  2(b)  quali¬ 
tatively  supports  this  expectation  over  the  displayed  range 
of  n. 

As  shown  in  Fig.  2(c),  the  ratio  of  molecular  and 
atomic  line  strengths  ranges  from  «1%  —  5%.  By  taking 
our  peak  atomic  density  into  consideration,  this  agrees 
quite  well  with  the  relative  signal  strengths  found  in 
Refs.  [2,6].  One  may  expect  the  molecular-signal  ratio  to 
scale  with  the  probability  of  finding  a  5Si/2  atom  within 
the  Rydberg-atom  volume  (which  scales  as  n 6),  corre¬ 
sponding  to  an  increase  of  about  a  factor  of  2.5  over  the 
n  range  in  Fig.  2(c).  The  observed  increase  is  about  a 
factor  of  5.  This  enhancement  is  likely  due  to  a  Rydberg 
excitation  blockade  caused  by  electrostatic  Rydberg-atom 
interactions  [17,18],  which  suppresses  the  atomic  line 
[19].  Since  the  blockade’s  effectiveness  increases  with  n , 
the  molecular-signal  ratio  scales  faster  than  n6.  This 
interpretation  is  corroborated  by  the  atomic-signal 
strength,  which  drops  by  a  factor  of  4  over  the  n  range 
in  Fig.  2(c).  In  the  absence  of  an  excitation  blockade,  the 
atomic  signal  would  drop  as  n~ 3,  i.e.,  by  only  a  factor  of 
2  in  Fig.  2(c). 

Molecular  states  of  the  nDj  +  5Si/2  type  exhibit  a 
transition  between  Hund’s  case  (a)  at  low  n  and  Hund’s 
case  (c)  at  high  n .  Most  of  our  data  trend  towards  Hund’s 
case  (c),  where  the  dominant  molecular  potentials  carry  spin- 
dependent  factors  |  {nif  =  0,  ms  =  ±l/2|j,  rrij  =  ±1/2)  |2 
(detailed  model  below);  these  factors  are  1)  for 

j  =  €-  1/2  and  (f  +  1)/(2*?±1)  for  j  =  jf  +  1/2. 
Hence,  in  Hund’s  case  (c)  (high  n)  the  v  =  0  binding-energy 
ratio  for  nDj  ±  5 Sy2  j  =  3/2  and  5/2  molecules  should  be 
about  2/3.  For  decreasing  n,  the  fine  stmcture  splitting 
increases  as  n~ 3,  while  the  molecular  binding  energy 
increases  as  n~6.  The  system  then  trends  towards  Hund’s 
case  (a),  in  which  the  v  —  0  binding  energy  ratio  changes 
from  2/3  to  1  (see  below).  Therefore,  the  v  —  0  binding 
energy  ratio  is  a  convenient  experimental  measure  to 
characterize  the  system.  Figure  3  shows  spectra  of  the 
two  fine  stmcture  components  of  37 D.  Molecular  peaks 
are  present  for  both  j  —  3/2  and  5/2,  with  binding  energies 
of  24  ±  3  and  26  ±  3  MHz,  respectively,  corresponding  to  a 
ratio  of  0.92  ±  0.15.  Since  this  is  significantly  larger  than 
2/3,  for  n  —  37  the  system  is  in  the  transition  regime 
between  Hund’s  cases  (c)  and  (a). 

Next  we  compare  experimental  data  with  a  calculation. 
We  consider  s-wave  triplet  scattering  and  omit  singlet 
scattering  [1].  For  low  electron  momenta  k ,  the  s-wave 
triplet  scattering  length  is  a]  (k)  «  as0  ±  |  ak,  where 
as0  —  aI(k  =  0)  is  the  zero-energy  scattering  length  and 
a  the  Rb  5Si/2  polarizability  [13].  We  fit  as0  to  match  the 
experimental  data.  In  calculations  not  presented  here  we 


FIG.  3.  Spectra  centered  on  the  37 Dj  atomic  Rydberg  lines  for 
j  —  3/2  (left)  and  j  =  5/2  (right;  the  same  as  in  Fig.  2). 
87Rb(37DJ  +55'1/2)(^  =  0)  molecular  signals  are  indicated  by 
vertical  dashed  lines  and  squares. 

have  verified  that  p- wave  scattering  strongly  influences  the 
inner  part  of  the  molecular  potentials  [20]  but  leaves  the 
outermost  well,  where  the  bound  states  relevant  to  our 
work  are  found,  largely  unaffected.  The  Rydberg  atom’s 
fine  stmcture  and  the  perturber’s  hyperfine  stmcture  are 
included,  because  they  are  on  the  same  order  or  larger  than 
the  molecular  binding.  The  much  smaller  hyperfine  stmc¬ 
ture  of  the  Rydberg  atom  is  ignored.  We  include  configu¬ 
ration  mixing  with  neighboring  Rydberg  manifolds  [5].  For 
the  5Si/2  atom  located  at  R  =  Zez,  the  Hamiltonian  is 

Hq  ±  2naTs  [fc(r)]<53(?  —  Zez)  •  S2  ±  — ^  ±  AS2  •  I2, 

(1) 

where  the  unperturbed  Hamiltonian  H0  includes  Rydberg 
quantum  defects  and  fine  stmcture  [16].  The  operators  Si 
and  S2  are  the  spins  of  the  Rydberg  electron  and  5 Sl/2 
atom,  respectively.  The  87Rb  5Si/2  atom  has  jf2  =  0,  a 
nuclear  spin  I2  with  /2  =  3/2,  and  a  hyperfine  parameter 
A  =  ft  x  3.4  GHz  (in  SI  units).  The  projector  Si  •  S2  ±  | 
has  the  eigenvalue  one  (zero)  for  the  triplet  (singlet)  states 
of  Si  and  S2,  enabling  only  triplet  scattering.  In  the 
classically  allowed  range  of  the  Rydberg  electron 
k  =  ^-l/(neff)1fteff)2)  ±  2/r  (atomic  units),  and  k  —  0 
elsewhere.  There,  neff)1  and  neff  2  are  the  effective  quantum 
numbers  of  the  Rydberg  states  coupled  by  the  scattering 
term.  Since  only  states  with  =  0  are  nonvanishing 
on  the  intemuclear  axis,  the  relevant  Hilbert  space  is 
limited  to  {\n,  jx,  m;1  =  ±1/2,  ms2  =  ±1/2,  mi2  = 
±1/2,  ±3/2)}.  Since  the  Hamiltonian  in  Eq.  (1)  conserves 
mk  :=  rriji  ±  ms2  ±  mi2,  the  space  breaks  up  into  separated 
subspaces  with  =  ±5/2,  ±3/2,  ±1/2,  where  the 
scattering  term  couples  states  with  the  same  m;1  ±  ms2, 
while  the  hyperfine  term  couples  states  with  the 
same  ms2  ±  mi2. 

The  Hamiltonian  in  Eq.  (1)  is  diagonalized,  resulting  in 
adiabatic  potential  surfaces  Vadi(Z)  and  electric  [<7z-(Z)] 
and  magnetic  [/l  (Z)]  dipole  moments  of  the  adiabatic 
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FIG.  4  (color  online).  (a),(b)  Deep  (solid  lines)  and  shallow 
(dashed  lines)  potentials  for  35 D3/2  and  35T>5/2-type  molecules, 
for  as o  =  -14a0,  and  vibrational  wave  functions  for  v  —  0,  1  in 
the  deep  potentials.  (c),(d)  Energy  levels  for  v  =  0,  1  in  the  deep 
(triangles)  and  shallow  (circles)  potentials  vs  n. 


FIG.  5  (color  online),  (a)  v  —  0  binding  energies  in  the  deep 
molecular  potentials  for  nD5/2  (filled  triangles),  nDy2  (open 
triangles),  and  the  D  fine  structure  splitting  (dashed  line)  vs  the 
effective  quantum  number.  Solid  lines  are  fits.  The  D3/2  energies 
are  fit  well  by  84.2  GHz /n^}3.  The  D5  j2  energies  do  not  exhibit  a 
global  scaling;  at  low  n  they  tend  to  scale  as  23  MHz/n^. 
(b)  Electric  dipole  moments  for  v  —  0  vs  n  for  the  deep  (triangles) 
and  shallow  (circles)  Vad(Z)  for  j  —  3/2  (open)  and  j  —  5/2 
(filled),  (c)  Magnetic  dipole  moments  for  the  same  states  as  in  (b). 


states  \i  is  an  arbitrary  label  for  the  Vad(Z)].  The  vibrational 
states  W i  v  and  their  wave  functions  ^>(Z)  are  found  by 
solving  the  Schrodinger  equation  with  potential  Vadi{Z) 
and  reduced  mass  87  amu/2,  and  their  electric  and  mag¬ 
netic  dipole  moments  are  di  v  =  f  l^i  l/(Z)j2di(Z)dZ  and 
Hi.v  =  J’|\I'I>(Z)|2//j(Z)</Z,  respectively. 

In  Figs.  4(a)  and  4(b),  we  show  all  potentials  Vad(Z)  for 
35 D5/2  and  35 D3/2  that  connect  to  the  F  =  2  hyperfine 
level  of  the  55^2  atom,  as  well  as  the  vibrational  states 
v  —  0,  1  of  the  deeper  potentials.  The  deep  and  shallow 
potentials  and  their  states  have  degeneracies  of  6  and  4, 
respectively.  The  shallow  potentials  are  due  to  hyperfine- 
induced  mixing  between  the  electronic  singlet  and  a  triplet 
state.  In  Figs.  4(c)  and  4(d),  we  show  the  vibrational 
energies  W i  v  for  v  =  0,  1  over  a  range  of  n  for  nD5 /2  + 
55^2  and  nD3/2  +  5S\/2  molecules. 

In  the  experiment,  we  excite  molecules  below  the  F  =  2 
asymptote,  for  which  the  states  in  the  deep  potential 
Vad{Z )  have  the  larger  degeneracy  and  are  easier  to  observe 
due  to  their  larger  binding  energies.  Hence,  we  adjust  as0 
so  that  the  v  —  0  levels  in  the  deep  V ad(Z )  match  the 
experimental  data  in  Figs.  1-3  and  find  as0  =  -14a0  d= 
0.5<2o  [see  Fig.  2(b)].  (For  37D5/2  the  v  =  0  binding  energy 
increases  by  about  4  MHz  when  changing  as0  from  -13.5 
to  —  14.5fl0*)  This  as0  value  lies  within  the  range  of 
published  values  -13  to  -19.48a0  [2,6,21-23].  The 
experimental  signals  at  about  half  the  v  =  0  binding 
energies,  seen  in  some  of  the  plots  in  Fig.  2(a),  may  cor¬ 
respond  to  levels  at  about  half  the  v  =  0  binding  energies 
in  Figs.  4(c)  and  4(d). 

In  the  high -n  limit  in  Fig.  5(a),  the  v  =  0  binding-energy 
ratio  for  D3/2  and  D5/2  approaches  2/3,  as  expected  for 
Hund’s  case  (c),  and  both  sets  of  binding  energies  approx¬ 
imately  scale  inversely  with  the  atomic  volume  (i.e.,  as 


At  low  n,  the  system  transitions  into  Hund’s  case  (a), 
where  the  binding-energy  ratio  changes  from  2/3  to  »  1. 
Also,  the  v  —  0  binding  energies  for  the  lower  (D3/2)  fine 
structure  level  exceed  the  fine  structure  coupling  and  keep 
scaling  as  while  those  for  the  D5/2  level  approach  the 
fine  structure  splitting  and  its  scaling  The  exper¬ 

imental  data  in  Fig.  3  are  in  the  transition  regime  between 
the  two  Hund’s  cases. 

For  nD  Rydberg  molecules  the  electric  dipole  moments 
dv= o  are  somewhat  smaller  than  for  nS  ones  [see  Fig.  5(b) 
and  Ref.  [7]].  As  a  result  of  the  transition  between  two 
Hund’s  cases,  the  dv=0  of  the  nD5/2- type  molecules  do  not 
exhibit  a  clear  scaling.  Those  for  the  nD3j2- type  molecules 
scale  close  to  n~^,  similar  to  S-type  ones  [7].  Noting  that 
dv= o  <  0.25 ea0  and  that  the  stray  electric  field  in  our 
experiment  is  <200  mV /  cm,  the  experimental  permanent- 
electric-dipole  shift  is  well  below  1  MHz.  The  magnetic 
moments  [Fig.  5(c)]  exhibit  a  gradual  change  that  reflects 
the  change  in  angular-momentum  coupling  behavior 
between  the  two  Hund’s  cases.  In  our  <  1  G  magnetic 
field  in  the  atom  trap,  the  molecular  lines  could  be 
broadened  over  a  range  <4  MHz,  while  the  atomic  lines 
have  \rrij\  <5/2  and  could  be  broadened  over  a  range 
<8  MHz.  These  estimates  agree  well  with  the  linewidths 
seen  in  the  spectra. 

In  summary,  we  have  observed  polar  Rydberg  molecules 
of  the  type  nDj  -f  55j/2-  The  binding  energies  are  larger 
than  those  of  previously  observed  nSi/2  +  5Si/2  molecules. 
The  molecules  undergo  a  transition  between  Hund’s  case 
(a)  at  low  n  to  Hund’s  case  (c)  at  high  n.  With  improved 
spectroscopic  resolution,  one  should  be  able  to  measure 
the  electric  dipole  moments  and  to  assign  higher-lying 
vibrational  levels  in  the  deep  potentials  in  Fig.  4  as  well  as 
levels  in  the  shallow  potentials.  Recently,  we  became 
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aware  of  related  work  on  nD  molecular  states  in  a  higher- 
magnetic-field  regime  [24]. 
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Interactions  between  the  excitations  cause  excitation  blockades,  short-  and  medium-range  order,  as  well 
as  dynamics  due  to  van-der-Waals  interaction.  All  these  topics  have  been  studied.  Two  papers  were 
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trapped,  for  the  first  time,  cold  circular  Rydberg  atoms  in  these  traps.  In  addition  to  being  able  to  study  the 
dynamics  of  these  highly  dipolar  atoms  in  the  magnetic  trap,  we  also  observed  a  collision-induced  loss  we 
have  traced  back  to  collisions  between  the  long-lived  circular  atoms  and  cold  ground  state  atoms.  These 
types  of  collisions  are  unexplored,  to  date.  In  this  setup  we  also  have  prepared  and  studied  D-type 
molecules  between  Rydberg  and  ground-state  atoms.  These  are  due  to  an  attractive  triplet  s-wave  low- 
energy  scattering  interaction  between  the  Rydberg  electron  and  a  ground-state  atom.  We  have  published 
two  papers,  one  in  PRL,  on  in  PRA. 
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